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We review how the renormalized force correlator A{u), the function computed in the functional RG field 
theory, can be measured directly in numerics and experiments on the dynamics of elastic manifolds in presence 
of pinning disorder. We show how this function can be computed analytically for a particle dragged through 
a 1-dimensional random-force landscape. The limit of small velocity allows to access the critical behavior at 
the depinning transition. For uncorrelated forces one finds three universality classes, corresponding to the three 
extreme value statistics, Gumbel, WeibuU, and Frechet. For each class we obtain analytically the universal func- 
tion A{u), the corrections to the critical force, and the joint probability distribution of avalanche sizes s and 
waiting times w. We find P{s) — P{w) for all three cases. All results are checked numerically. For a Brownian 
force landscape, known as the ABBM model, avalanche distributions and A(it) can be computed for any ve- 
locity. For 2-dimensional disorder, we perform large-scale numerical simulations to calculate the renormalized 
force correlator tensor Aij{u), and to extract the anisotropic scaling exponents (x > Ca- We also show how 
the Middleton theorem is violated. Our results are relevant for the record statistics of random sequences with 
linear trends, as encountered e.g. in some models of global warming. We give the joint distribution of the time 
s between two successive records and their difference in value w. 



I. INTRODUCTION 

Elastic objects driven through a disordered environ- 
ment are ubiquitous in nature, including magnets [1, 78], 
superconductors [2, 3], density waves [4, 5], wetting [6, 7], dry 
friction [8], dislocation [9], crack propagation [10], and earth- 
quake dynamics [11]. These phenomena can be studied by 
different theoretical approaches, including phenomenological 
arguments [2], mean field models [12], functional renormali- 
sation group for statics [13-39] and driven dynamics [40-47]. 
They were also studied with numerical techniques [48-51]. 
In several cases the experimental results seem to be in reason- 
able agreement with the theory (see [52] for vortex lattices. 
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FIG. 1: Dynamical shocks (avalanches): position of a particle Uw 
pulled by a spring, of varying equilibrium position w, in a one di- 
mensional random force landscape (with forces uniformly distributed 
between and 1). The quasi-static motion shows a succession of 
jumps, also called shocks. Decreasing the spring constant (the mass) 
from = 0.01 (red) over = 0.03 (green) to = 0.001 
(blue), the shocks become larger and larger. 



[53] for ferroelectric s , and [1] for magnetic interfaces), but 
some discrepancies are still manifest, at least with the sim- 
plest theories, in some cases, e.g. the depinning of the contact 
line of a fluid [6,7,47]. 

Recent theoretical progress allows not only for qualitative, 
but also for quantitative tests. On one hand, for interfaces, 
powerful algorithms now allow to find the exact depinning 
threshold and critical configuration on a cylinder [51, 54] and 
to study creep dynamics [55]. On the other hand the func- 
tional RG (FRG) has been extended beyond the lowest order 
(one loop), and it was shown that differences between statics 
and depinning become manifest only at two loops [13, 40], 




FIG. 2: (a): The measured correlator A(ii) for a particle pulled in 
a random potential (i.e. RB disorder) distributed uniformly in [0, 1], 
and rescaled such that A(0) = 1 and du\A{u)\ — 1. From 
bottom (which has Jg°° duA{u) ~ 0) to top the mass decreases from 
= 1 (red) to — 0.5 (green) to = 0.003 (blue), where 
it has (up to small corrections) converged to the fixed point A*(u). 
This demonstrates in d = the expected crossover from random- 
bond to random-field disorder (see text). Note that the fixed point 
has a cusp singularity at it = 0. 
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i.e. to second order in an expansion in rf = 4 — e where d is 
the internal dimension of the manifold. Such differences ap- 
pear for instance in the roughness exponent (. The FRG is 
a field theoretic tool for disordered systems, which captures 
the complex glassy physics of numerous metastable states at 
the expense of introducing, rather than a single coupling as 
in standard critical phenomena, a function, A(u), of the dis- 
placement field u, which flows to a fixed point (FP) A*{u). 
This FP is non-analytic, as is the effective action of the theory. 
The non-analyticity is a rather unconventional feature, and the 
validity of the approach has been questioned: one could argue 
that, although A*(ii) is perturbative near d = 4, the second 
derivative A"(0) has gone to infinity, hence we have left the 
domain controlled by perturbation theory. To put this criti- 
cism to rest, one first shows how an observable of the exper- 
imental system can be defined, which is identical to the field 
theoretic disorder correlator A{u): the idea [56] is to add a 
quadratic confining potential to the system, which formally 
acts as a mass for the elastic modes of the interface. The dis- 
order correlator A(w,), can then be measured directly as the 
second moment of the interface displacement. 

This method has been used in a numerical simulation of in- 
terfaces in a disordered magnet, to compute numerically the 
zero-temperature FRG fixed-point function A(u) in the stat- 
ics, for interfaces (N = 1) using powerful exact minimization 
algorithms [59]. A variety of disorder types, random-bond, 
random-field and periodic disorder were studied in various di- 
mensions = 0, 1, 2, 3. The results are close to 1-loop predic- 
tions and deviations are consistent with 2-loop FRG. The most 
important feature, namely a linear cusp in A(u) was clearly 
seen. These results come in strong support for the underly- 
ing hypothesis of a non-analytic field theory, perturbatively 
accessible in a d = 4 — e expansion. 

In our Letter [60], we have extended the method of [56, 59] 
to driven systems. In particular it allows to measure the FRG 
fixed-point function A* (u) near the depinning transition at ve- 
locity w = 0+. The form of this fixed point was obtained to 
one loop [42, 43], but the remarkable fact is that it is only 
to two loop order that it differs from the static fixed point 
[13, 40]. Shortly after, this tiny difference was measured, be- 
yond statistical uncertainties, together with the predicted lin- 
ear cusp, in a numerical study [61] of a line driven in a one 
dimensional medium. 

In this situation, it is useful to find a simple model, which 
can be solved analytically, and exhibits many features of the 
more complicated situation. Such a model is a particle in a 
random force landscape, pulled by a moving spring. Since 
there is no internal degree of freedom, it is the d ^ limit 
of the depinning fixed point for interfaces. Similar d = 
toy models were very useful in the study of the statics. The 
universality classes there were found to be parameterized by 
the exponent ( > 1, with a (presumably unique) univer- 
sal fixed-point function A{u) = —R"{u) in each case, with 
R{u) ^ at large \u\, with 7 = 4 — 4/^. Only in some 
cases this function was obtained analytically, e.g. for the case 
C = 4/3 of the Sinai model [56], which corresponds to the 
random-field disorder class. It is thus quite interesting to ob- 
tain the corresponding results for the depinning fixed point. 



Of course another application of the d = model is driving 
a fixed-size manifold over very large distance, it eventually 
behaves again as a particle with some effective random force 
landscape. 

The aim of this paper is to give a detailed account of the 
results summarized in [60], and to present some new ones. 
We summarize the basic ideas in section II, relegating details 
of the field-theoretic derivation in appendix A. 

In section III, we give a detailed derivation of the analytical 
results for the d = particle model, focussing there on uncor- 
related random forces. This yields an exponent C in a continu- 
ous range < ( < 2 with corresponding fixed-point functions 
A(u) easier to compute than in the statics. The result depends 
on the tail of the distribution of the local random force, i.e. 
we find the three main universality classes of extreme value 
statistics: The Gumbel, WeibuU and Frechet distribution. 

In section IV, we calculate analytically the joint avalanche 
size s and waiting-time w distribution. The avalanche distri- 
bution was computed recently [62] in a d — 4 — e expan- 
sion using FRG, and the present results hence correspond to 
the d = Umit. There it was shown how the avalanche- 
distribution is related to the non-analyticity of the set of aU 
cumulants of the displacement field. 

Interestingly, the problem of a particle driven in one di- 
mension is related to the so-called record statistics [73, 74]. 
If the particle is pulled by a spring the problem is related to 
record statistics for random sequences with linear trends [76], 
whose interest has been revived in the context of global warm- 
ing models [75]. Translated into the language of records, we 
obtain in Section IV, the joint distribution of the time s be- 
tween two successive records and their difference in value, w, 
for a sequence of i.i.d. random variables with trends, i.e. vari- 
ables Yn = Xn + cn where X„ are i.i.d. random variables, 
and c a drift. These results are obtained for the three classes 
of extremal statistics. 

In section V, we check some of our above results numeri- 
cally; we also study numerically a particle driven at non-zero 
velocity v > 0, and find that the velocity smoothens the cusp 
in the force correlator. 

In section VI we consider long-range correlated random 
force landscapes, specifically the case of a Brownian force. 
For this model, known as the ABBM model for domain wall 
motion, remarkably, the stationary distribution of instanta- 
neous velocities can be computed [77] for any non-zero aver- 
age driving velocity v > 0. From that we obtain A(u) for any 
V > 0. We also compute the quasi-static A(u) and avalanche 
distribution. It matches with the limit v 0+, and shows 
how the cusp is smoothened at t; > 0. 

In section VII we summarize some known results and some 
new ones, common to record statistics and to the present 
model of a driven particle, either with no mass (fixed-force 
driving, symmetric records), or with a mass (fixed-velocity 
driving, records with drifts). In particular we study in detail 
the record statistics for a Levy-walk landscape with drift. 

Finally we address the outstanding question of the depin- 
ning for systems which can move in more than one direction, 
also termed "N > 1", with N the number of components 
of the displacement field u (e.g. = 2 for a line moving 
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in three dimension). In particular there is still no satisfac- 
tory field theoretic description for this case based on FRG. 
This question was studied by Ertag and Kardar [63], but they 
made the approximation that the disorder correlator only de- 
pends on the direction in which the system is driven. Consid- 
ering two manifolds which are driven on trajectories far apart 
in the transversal direction, their renormalized disorder corre- 
lators should be independent, questioning the assumptions in 
the Erta§-Kardar approach. We have studied this situation in 
the field theory [66] (see also some study at > in [67]), 
but consistent and stable solution of the fixed-point equations 
seem quite complicated and are still lacking. In this situa- 
tion it is important to have some numerical results as guide 
for the analytical treatment. In section VIII, we therefore dis- 
cuss the changes necessary to study an elastic manifold driven 
through a higher-dimensional random environment, and com- 
plement this in section IX by a numerical study of a particle 
dragged through a random energy landscape. Especially, we 
show numerically, that the scahng exponents ("roughness") in 
the direction of the driving (x and perpendicular to it (y are 
different, and satisfy Q > Cy We also find, that the cross- 
correlator in the transversal direction (i.e. the force correla- 
tor between forces in the direction x of the driving, and its 
transversal one y, measured as a function of the transversal 
distance Uy, Axy{ux = 0, %), is non- vanishing. 

II. SUMMARY OF THE METHOD 

A. General framework 

We consider the equation of motion for the over-damped 
dynamics of an elastic manifold parameterized by its time- 
dependent displacement field u{x, t): 

r]dtu{x,t) = Fx[u{t);w{t)\ (1) 
Fx[u;w] = m^{w - u{x)) + cV^u(a;) -|- F{x,u{x)) 

where Fx[u{t)]w{t)] is the total force exerted on the mani- 
fold (we note u{t) = {u{x, t)}x^^d the manifold configura- 
tion, X being its d-dimensional internal coordinate); r] is the 
friction coefficient and c the elastic constant. Here at the bare 
level, the random pinning force \sF{x,u) = —duV{x,u) and 
the random potential V has correlations y (0, x)V{u, x') = 
Ro{u)S'^(x — x'). We consider first random-bond bare disor- 
der with a short-ranged Ro{u). We have added a harmonic 
coupling to an external variable w{t), a given function of time 
(in most cases we choose it uniform in x). This is the sim- 
plest generalization of the statics, where w{t) = w is time- 
independent. It is useful to define the fixed-w energy 

Hy^iu] = / d'^x —{uix) - u>)2 + V{x, u{x)) (2) 

associated to the force Fx[u; w] = —^^j^f^- If w{t) is an in- 
creasing function of t the model represents an elastic manifold 
"pulled" by a spring. 

We first describe qualitatively how to measure the FRG 
functions and later justify why the relation is expected to 



be exact. We are interested in the observable w{t) — {u{t)) 
where u{t) = L~'^ J d'^xu{x, t) is the center of mass po- 
sition, and (...) denotes thermal averages, i.e. the ground 
state at zero temperature. It represents the shift between the 
translationally averaged displacement and the center of the 
well, i.e. the extension of the spring. It is thus proportional 
to the pulling force on the manifold, hence to the transla- 
tionally averaged pinning force minus the friction force, i.e. 
w{t) — u{t) = rnr'^{riv{t) — j,^ F(x, u(x, t))) (if we use peri- 
odic boundary conditions inside the manifold). Of particular 
interest are: 

wit) - {u{t)) = m-\fay{t) (3) 
{w{t) - {u{t))){w{t') ^ {u{t'}) f = m-^L-'^DUt,t') 

where cormected means w.r.t. the double average (...). If we 
consider a function w{t) such that dw{t) /dt > 0, one can also 
write: 

D^{t,t') = Ay,{w{t),w{t')) (4) 

As written, the function A^, may in general depend on the 
history w{t). However we expect that for fixed L,m and 
slow enough w{t), e.g. w{t) = vt with v 0"*", one 
has A„(w(t),w;(f')) A{w{t) - w{t')). This function 
A(ii; — w'), which is independent of the process w{t), is the 
one defined in the F.T.. The derivation of this property is given 
in Appendix A to which we refer the reader for technical de- 
tails. Note that we are discussing now = 1 systems (in- 
terfaces), subtleties related to > 1 are discussed in section 
VIII. 

Let us now describe T = depinning. Quasi-static de- 
pinning is studied as the limiting case where dw/dt 0+. 
The quasi-static motion can be described as follows (in the 
continuum model). One starts in a metastable state uo{x) for 
a given w = wq, i.e. a zero-force state Fj;{uo{x);w) = 
which is a local minimum of Hj^,„ [u] with a positive barrier. 
One then increases w. For smooth short-scale disorder, the 
resulting deformation of u{x) is smooth. At some w = wi, 
the barrier vanishes. For w = wf the manifold moves down- 
ward in energy until it is blocked again in a metastable state 
fii(x) which again is a local minimum of We are 

interested in the center of mass (i.e. translationally averaged) 
displacement u = L~'^ J d'^xu{x). The above process de- 
fines a function u{w) which exhibits jumps at the set Wi. Note 
that time has disappeared: evolution is only used to find the 
next location. The first two cumulants 

w — u{w) = m~^fc (5) 
(w - u{w)){w' - u{w'f = m~'^L^'^A{w - w') (6) 

allow a direct determination of the averaged (m-dependent) 
critical force fc and of A{w). Note that u{w) a priori de- 
pends on the initial condition and on its orbit but at fixed m 
one expects an averaging effect when w is moved over a large 
region. This is further discussed below. Note that the defi- 
nition of the (finite size) critical force is very delicate in the 
thermodynamic limit [38]. 
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Elastic systems driven by a spring and stick-slip type mo- 
tion were studied before, e.g. in the context of dry friction. 
The force fluctuations, and jump distribution were studied nu- 
merically for a string driven in a random potential [64]. How- 
ever, the precise connection to quantities defined and com- 
puted in the field theory has to our knowledge not been made. 
The dependence in m for small m predicted by FRG, A{w) = 
rrf^'^^K{'wmr^) is consistent with observations of [64] but 
the resulting A (it;) has never been measured. Fully con- 
nected mean-field models of depinning also reduce to a par- 
ticle pulled by a spring, together with some self-consistency 
condition. Ref. [43] discusses related issues in an expansion 
around mean field. As discussed below, our main remarks 
here are much more general, independent of any approxima- 
tion scheme, and provide a rather simple and transparent way 
to attack the problem. 

Note that the manifold in the harmonic well can be ap- 
proximated by {L/L„iY roughly independent pieces with 
Lm ~ 1/m. The motion of each piece over large distances 
resembles the one of a particle, i.e. a d = model, but with a 
rescaled unit of distance in the u direction, Um ^ L^n ^ rn^^. 
The "effective force" landscape seen by each piece becomes 
uncorrelated on such distances, and its amplitude scales as 
Fm ~ rn^Um- Hence one is in a bulk regime not dominated 
by extremes, i.e. A{w) probes only motion over order one 
unit. It is easy to check that an arbitrary initial condition joins 
the common unique orbit after about one correlation length. 
Hence the d = model suggests that starting the quasi-static 
motion in uq and driving the manifold over w ^ L should 
then result in all orbits either converging or having statistically 
identical properties. Note that if the manifold is driven over 
more than L^, a crossover to d = behavior and extremal 
statistics occurs, as studied in the next Section. 

The averaged critical force, defined in (5), should, for d > 
0, go to a finite limit, with fc{m) — + Bm?^^ from finite 
size scaling. This has been recently tested in the numerics 
[61]. Although fc is not universal and depends on short-scale 
details, one easily sees that —mdmfci'm) depends only on one 
unknown scale. We note that the definition (5) coincides with 
the one proposed recently as the maximum depinning force 
for all configurations having the same center of mass uq [38]. 
Since u—w is a fluctuating variable of order {L j Ljn)~'^l'^ , the 
two definitions should coincide in the limit where L ^ oo, 
before m — > 0. The one point distribution of the critical force 
is obtained from the distribution of w -~ u{'w), and to one loop 
is identical to the one obtained in [38] provided one uses there 
the massive scheme. 

Let us now recall the field-theory predictions: The FRG 
equation at 2-loop order for the (rescaled) force correlator are 
[13]: 



A{u) 




FIG. 3: Fixed-point functions A*{u) for random-bond (RB) and 
random-field (RF) disorder (arbitrary scale). 



FIG. 4: Construction of u(w) in d = 0, for the pinning force F(u) 
(bold black line). The two quasi-static motions driven to the right and 
to the left are indicated by red and green arrows, and exhibit jumps 
("dynamical shocks"). The position of the shocks in the statics is 
shown for comparison, based on the Maxwell construction (equiva- 
lence of light blue and yellow areas, both bright in black and white). 
The critical force is 1/ (2A/) times the area bounded by the hull of 
the construction. 

A = — 1 describes the statics, and A = 1 the depinning. For 
the statics, A(u) admits a potential solution (random-bond 
universality class) A(it) — —R"{u), with R{u) decaying to 
as It oo, as illustrated on figure 3. This implies that the in- 
tegral du A(u) remains unrenormalized. However at de- 
pinning, it flows, and no potential solution exists. For a large 
class of bare disorder, the model should renormalize, as m 
decreases, to the random-field fixed-point solution, A^p(m), 
which is monotonically decaying and strictly positive. For the 
d = toy-model discussed in the next section, this crossover 
is nicely seen in our numerical simulation with decreasing to, 
as is illustrated on figure 2. Therefore, in the following, we 
can focus on the random-field universality class, i.e. short- 
range correlated random forces. 



-mdmA{u) = (e - 20 A(w) + CuA'{u) HI. PARTICLE IN SHORT-RANGE RANDOM-FORCE 

2 ri ^ LANDSCAPE (d = 0): EXACT RESULTS 

d I 1 r . / N . ,„Ni2 1 



^ ' \2rA^,^ _ ArnM^\ AVn+^2A^,^\tlonof motion 



We now study the model in = 0, i.e. a particle with equa- 

+^ {A'(m)2[A(u) - A(0)]+AA'(0+)2A(w)}t' 

(7) 7jdtU = m^{w-u) + F{u) . (8) 
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FIG. 5: Construction of u{w) in d = for the forward motion in the 
discretized model. The vertical lines are the force barriers, the (red) 
increasing lines the spring force rn^(u — w). A particle moves from 
left to right, until it is stopped by a barrier (when the lines for spring 
force and barrier forces intersect). 



In the quasi-static limit where w is increased slower than 
any other time-scale in the problem, the zero force condition 
F{u) = m^(u — w) determines u{w) for each w, starting 
from some initial condition. The graphical construction of 
u{w) is well known from studies of dry friction [? ]. When 
there are several roots one must follow the root as indicated in 
Fig. 4, where F{u) is plotted versus m?{u — w). This results 
in jumps and different paths, v}(w) and u^(yS) respectively 
for motion to the right (forward) and to the left (backward). 
Let us call A the area of this hysteresis loop (the area of all 
colored/shaded regions in Fig. 4). It is the total work of the 
friction force when moving the center of the harmonic well 
quasi-statically once forth and back, i.e. the total dissipated 
energy. The above definition of the averaged critical force (5), 
assuming the landscape statistics to be translationally invari- 
ant and that one can replace disorder averages by translational 
ones over a large width M (which certainly holds if force cor- 
relations are short-range correlated), gives 



A. Short-range correlated force: a discrete model 

Let us now consider the asymptotic forward process 
u(w) := u^{w), defined in general as the smallest root of 
the equation F{u) — m?{u — w). The set of points (w, F{u)) 
is the lightened portion of the F{u) curve, the rest being the 
shadow. If one starts in the shadow at wo one joins at wi the 
asymptotic process for all w > wi. The difference wi — wq 
is finite for finite m hence we will only study the asymptotic 
process. Note that the area of the shadow per unit length is the 
critical force. 

We now study short-range correlated random force land- 
scapes. In the limit of interest, m — > 0, the scale of jumps 
becomes large, and the finite range should be unimportant. 
Hence it is equivalent, and more convenient, to consider a dis- 
crete model, u being integers. The variable w can be kept real. 
One considers a discrete landscape F{u) = Fi independently 
distributed with P{F), and i integer variable. The process 
u{w) is then defined on integers. Its definition is shown in 
Fig. 5. 

Let us compute for w' > w the following joint probabiU- 
ties: 

P^U, F) := Pmh{u{w) = j and F^ = F) 

(11) 

Pw-w'U, F;j', F') := Proh(u{w) - j and F, - F 



We define; 



and u{w') = j' and Fj> = F' (12) 



/'+00 />F 

H{F) = / P(/)d/ = 1 - / P(/)d/ (13) 

J F J-oo 



2 r^i 



fc = m (w - Uw) = / dw {w~ u„) (9) 



Since foTu{w) = j to hold, one must have all Fk > rn? [k- 
w) for all k < j and Fj < m^{j — w), see Fig. 5, one has; 



Hence, subtracting the two paths gives 



,2 fM 



fc — f}) = lim tttt / dw (u^(w) — uUw)) 



= lim — — 

A/^oo 2M 



(10) 



where we have used j u Aw — J wdu and < 0. One 
can check that for m ^ this definition of fc becomes 
identical to the one on a cylinder, fii, which for a particle 
(d = 0) is 2/rf ^ jI~ f^= max„ F{u) - min„ F{u) with 
2fdM — limm^o A{in). Since A depends on the starting 
point, this definition holds after a complete tour, where the 
maximum (minimal) pinning force was selected. One can also 
compare with the definition of shocks in the statics. There, the 
effective potential is a continuous function of w. Therefore, 
when making a jump, the integral over the force must be zero, 
which amounts to the Maxwell construction of figure 4. 



P^{j,F) = P{F)e{m^{j -w) - F) W Him^ik-w)) 



fc=-c 



i-1 



P^ij) = {l-Him^{j-w))) Y[ H{m\k~w)) A^^) 



k— — oo 

where the first line integrated over F yields the second (and 
6{x) denotes the unit step function). One easily checks that 

Next for / > j one has; 

P„;^,(j, = P{F)e{m\j -w)^F) 

x9{F~m^{j^w')) Yl H{m^{k~w)) (15) 



fc — — OO 



'.P{F')d{m'^{j' -w')-F') Y[ H(m^{k-w')) 

fc=i+i 
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with the convention that for j' 

=3 



j + 1, the factor define: 



nCZ+i H{w?{k - w')) = 1, and for j = j': 



Pw;n,'{j, F,j, F') = S{F - F')P{F)e{m\3 - w') - F) 
X H{m^{k-w)) (16) 

k=—oo 

Integrating over the forces we obtain the 1 -point and 2-point 
probability for the process u{w). They read, in compact no- 
tations with iJ^ := H{w?{k — w)), and non-zero only for 
/ > j: 



P^{3) = {l-HJ) n HI 

k=—oo 



(17) 



,7-1 j'-l 

p^.,^,{j,f)={Hf -Hyxi-Hp') n n 

k= — oo k=j+l 

+Sjj,{l-Hf) n (18) 



k— — oo 



l-HV 



/-I 



k=j+l 



(19) 



with the convention that for / = j, the factor ni=/+i ^k' — 
0. Using that E/>, (1 - ^/') ni=/+i = 1 one checks 
the normalization Ej/>j Pw;w'{j,j') = Pw{j)- 



B. Continuum limit for 1-point distribution 

For small m the continuum limit can be taken: 

P^(j)« [l-H{m\j-w)) 

xexp(^J dy In H{m'^ {y - w)j (20) 

/m'^{j-w) 
d/P(/) 
-OO 



: exp M d?/ ln(l - J 



m (y—w) 



m? {j — w) 



d/P(/) exp 



^ — OO J — 



P{l)dl) 

(y—w) 



P{l)dl . 



The last step is justified if the result is indeed dominated by the 
tail of P{f) for / negative, which is at the heart of extremal 
statistics. That this is indeed true is justified a posteriori. 

The quantity [w — u{w)] = {w—j) is the "local", i.e. 
fluctuating critical force, and its disorder average is fc{m) = 
Tn^[w — u{w)]. Its distribution can be obtained from the 1- 
point distribution Pw{j)- To rewrite (21) in a simpler form we 



P(/)d/ (21) 

-OO 

a,„(— oo) := . (22) 

Note that (+oo) = -|-oo. The 1-point distribution can thus 
be rewritten as: 

J'«,(i)di =e-'^-«da^(i). (23) 
Hence one can rewrite the first moment in the form: 



w — u{w) 



dj Pw{j){w- j) 



-OO 

OO 



dja;(j)e-»»W(w-j) (24) 



This means that the quantity a has a simple exponential dis- 
tribution, hence one needs to invert the relation, i.e. find j as 
a function of a 



j =j{a;w) 



a-wij) = a , 



(25) 



and use that a has an exponential distribution, to get any aver- 
age, e.g. 



{w - u{w))P = / da e-^iw - j(a; w)f (26) 
Jo 



for any p. 

C. Distribution of critical force : the different disorder classes 

We now obtain the universaUty classes for the 1-point dis- 
tribution of the process u{w), i.e. for the distribution of critical 
forces. We define b, a rescaled version of a: 

b{m^{j - w)) := n?a (27) 

Using (21), it can be written as 

h{x) = e-'=*(-) = r dy (28) 

J —OO J —OO 

The condition defines (j)- We thus need to invert this rela- 
tion to get j as a function of a. We do this for the three main 
disorder classes below. 



1. Gumbel class (class la) 

The first class contains distributions P{F) with unbounded 
support and decaying exponentially fast at F ^ — oo (in some 
broad sense defined below). One then finds the Gumbel distri- 
bution for the critical force, hence we call this class the Gum- 
bel class. 
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Let us invert relation (27) and assume that the following 
expansion holds at small m: 



m^{j — w) = (3 ^(In — Ina) 

In a 



-fcim) - -57 



/3'(-/0(m)) 



+ (r^)"(ln^) 



1 (lna)2 
2) 2 



/0(m) :=-6-i(m2)^-r^(ln^) 



(29) 
(30) 



where we know that a is a fluctuating number of order one 
with an exponential distribution, P{a) = e~°'9{a). That gives 
the distribution of the variable (j — w). In particular since 



^0 



da e"" In a = 7b = 0.577216 . . . , 



(31) 



we obtain the asymptotic s of the (averaged) critical force as: 



/e(m) = f,{m) 



/3'(-/°M) 



(32) 



For this asymptotics to hold, the ratio of successive terms has 
to go to zero, equivalent to 



lim —ln\{p-^y{z)\ . 

z^oo az 



(33) 



This defines the Gumbel class la, together with the fact 
that the support is unbounded. An example is lnP(/) = 
as / ^ -00. Then: 



0{x) w^^oo A{-xy 



7A 



(34) 
(35) 

(36) 



Hence class-la condition is satisfied for any 7, even for 7 < 1. 
Note that pre-exponential algebraic factors to not change the 
result. The critical force becomes 



m^{u{'w) — w) 
= rn? {j — w) 



1/7 



+ 



Ina 



In m~-^ \ 1- 



o 



(lna)2 



(lnm~^) 



(37) 



Defining the fluctuating critical force as fc{m) = 'm?{w 
u{w)) we find 



(38) 



where c = —Ina has a Gumbel distribution P(c) = 
e~'^ exp(— e~'^) on the real axis c e] — 00, oo[, with c = je 
and 



f°{m) = A"-(lnm"2)^ . 



(39) 



One also finds 

(/?-i)'(lnm-2) 



Pm = 



(40) 



Since as confirmed below, pm ^ "i"'' is the unique scale 
appearing also in the second cumulant (the disorder correlator 
defined in FRG), we can identify for class la: 



m ^ = 171 ^(Inm ^ . 



(41) 



Hence C = 2, with additional logarithmic corrections, i.e. ( = 
2+ for 7 < 1 and C = 2" for 7 > 1. 



2. Class lb: bounded support with exponential singularity 

An example of this class is: 

P{x) = A/{x + xoy (42) 

l3-\y) = -xo + {y/A)-^''< (43) 

h{x) ~ e-^/(^+=^'')^6l(x + ,xo) (44) 

P{f ) ~ e-^/(/+^°)^6'(/ + .To) (45) 

with 7 > 0. One sees that the condition (33) is obeyed. Hence 
this is still the Gumbel class, although we introduce a distinc- 
tion for convenience. A similar asymptotics can then be per- 
formed: 



"^^(i -w) = -Xo + 



A 



In) 



ln( 



7 (lnm-2)i+^ 



(46) 



Apart from the bound xq, the result is the same as in class I 
with 7 — > —7 (in the exponents only). Hence one finds again: 

fc{m) = f,{m) + cw?p^ (47) 

where c = —Ina has a Gumbel distribution P(c) = 
e~'^ exp(— e~°) on the real axis c e] — cjo, C!o[ with c = 7^ 
and: 



/°(m) = xo- A^'^{\iim-'^)-^ 



(48) 



One also finds the characteristic scale /9„ 



m 



l/(7A-i/^m2(lnm-2)^+^) , (49) 



hence C = 2+. 



3. Algebraic bounded support: Weibull class (class III) 

The Weibull class, or class III applies for a force distribu- 
tion with bounded support (from below) and algebraic behav- 
ior near the edge. An example is 



P{f) = Aif + xofe{f + xo) 
b{x) = A{x + xo)°'0{x + Xo) 

a ^ 2 + a 

A = A/[{2 + a){l + a)] . 



(50) 
(51) 
(52) 
(53) 
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One must have a > — 1 hence a > 1. The box distribution 
corresponds to a = 0, i.e. a = 2. Here 



(54) 



with y >0. Hence analogously to (29) and (30) we find from 
(27) 



m^(u(w) —w)= m^{j — w) = b ^{m^a) 
= -xo 



^) aV« , (55) 



where we recall that a is a random variable with distribu- 
tion P{a) = e~°'9{a). For a ^ oo one recovers the Gum- 
bel class. Hence one finds for the fluctuating critical force 

/c(m) = 'm?{u{w) — w) = rn?{j — w): 



fc{m) = fc{m) + cm'^pr, 



(56) 



where now c = — o^/" has a WeibuU distribution P(c) = 
a(— c)"-i exp(— (— c)") with parameter a on the negative 
real axis c g] — oo, 0[ and 

/°(m) = xo . (57) 

The averaged critical force is 

/e(m) = xo - mi r(l + i) + . . . , (58) 

where we have used that dae~°'a}/°' = T[l + i). One 
also finds that 



Pm 



^-l/a^-2(l-i) 



C = 2-- 
a 

with 1 < a < oo, hence < C < 2. 



(59) 
(60) 



4. Frechet class ( class II) 

The Frechet class, or class II, is relevant for force distri- 
butions with large fluctuations, i.e. algebraic tails on an un- 
bounded support. An example is P{f ) ~ A{—f)~°'Q{—f), 
a > 1. One has b{x) = A{—x)~°' with a = a — 2 and 

A = A{a - 2) (a - 1). Since (3{x) = - In A + a ln(-a;) and 



~ —6' 



one checks that the class-I condition (33) is 



not fulfilled. Let us first study a > 2, i.e. a > 0: 



(61) 



Hence one finds for the fluctuating critical force fdm) 
m^{u{w) — w) = m^{j — w): 



fc{m) = cm^pr, 



(62) 



where now c = a~^l°' has a Frechet distribution P(c) = 
ac-^-i exp(— c^") with parameter a > on the positive 
real axis c e]0, oo[ and the average critical force is: 



1 



/c(m) = r(l - -)m'p„ 



(63) 



where we have used c = da e "a ^Z" = r(l — and 
p„ = Ai/«m-2(i+4) . (64) 



This corresponds to a roughness exponent 



(65) 



Note that for a < 1 (a < 3, C > 4) the average critical force 
is infinite. 

We will see that the Frechet class is a bit pathologic in the 
sense that A(0) = oo for a < 2. More generally cumulants 
of order larger than a are infinite, i.e. they are associated to 
a probability distribution with fat tails. This implies as usual 
that these quantities are dominated by the largest events, hence 
they are sensitive to how the continuous limit is constructed 
from the discrete model. For a < 2 the integral in (28) is 
divergent at its lower bound, hence undefined without a cutoff. 



5. Comparison with extremal statistics 

Until now it seems that we have recovered the standard ex- 
tremal statistics classes for the distribution of the local critical 
force. On one hand this is not surprising, since one expects, 
qualitatively, that the critical force for two independent con- 
secutive regions in one dimension, be the maximum of the 
ones for each single region. This is certainly an exact state- 
ment for the zero-mass case recalled in Section VII. Here, it 
is quite consistent with the identification of ■m'^{u{w) — w) 
as a fluctuating threshold force, a fact which maybe was not 
obvious from the start. However, note that the value of the 
parameter a of the extremal statistics classes is shifted by one 
from the value it takes if one models the total critical force as 
the extremal one /c = minj({/^}) over ~ m~'^ indepen- 
dent regions, each with its critical force fl distributed with 
Pc{f)- In that case one has: 



Proba(/c > x) 



Pc{f) 



(66) 



which can be written equivalentiy as: 

J — oo 



(67) 



where a is a random variable of order one with an exponential 
distribution, P{a) = e~'^9{a). A comparison with (28) shows 
that the effective critical force in an independent region should 
be chosen with a distribution with a tail: 



=(/) = 



P{f)df 



(68) 
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for large negative /. Hence there is a first coarse graining 
which transforms the tail of P{f) into the tail of Pdf)- Then 
one can think of the resulting critical force as the maximum 
over N ~ independent random variables distributed 

with Pc{f)- Applied to WeibuU and Frechet classes, (68) in- 
deed accounts for the shift of the index a by one. 



D. 2-point probabilities and the FRG correlator A(w) 

In addition to recovering the three extremal statistics 
classes, which, as explained above, is not surprising if one 
thinks in terms of coarse grained independent random vari- 
ables, there is a second, more remarkable property. We find 
here that the 2-point correlation 



{w — u{w)){w' — u{w')) 



= J di dj'P^,^/(i, - j)iw' - j) -{w- j) {w' - f) 
= m-'^A{w - w') (69) 
takes, for all three classes, the following form at small m: 

A{w) = mV^A(«;/p„) , (70) 

where the scale Pm ~ to~^ is the one identified in each case in 
the previous section, and the fixed-point function A{w) only 
depends on the universahty class: it is unique and identical for 
all members of class I (la and lb), continuously depending on 
a (hence Q for classes II and III. 

Let us now give the joint probability in the continuum limit. 
From (18) and as was done to arrive to (23) one finds for w' > 
w: 

Pu,;u,'{j,f) = K{j) - a^,{j)] a'^,{j') (71) 

Let us check the normahzations. Using that 

Jp' dj'a'^,{j')e-''^''^^'^ = e-"»'(^) one obtains (writ- 
ing separately the two contributions): 



) 

= Pn,{i) (72) 
A similar trick yields 

/oo 
djPn.;^'{j,j')=P^'{j'). (73) 
-co 



Note also that 



(74) 



For all three classes one finds, starting from (37), (56) and 
(62) respectively 

j-w = -^^-pmc{a), (75) 



with either c(a) = — Ina (class I), c(a) = —0^/°" (class III) or 
c{a) = a~^/" (class II). Since the constant piece in (75) dis- 
appears when computing the connected moments, to compute 
A defined in (69) we can simply write 

j -w = -pmc{a) , f -w' = -pm.c{a') (76) 

and to obtain the rescaled function. A, we can further set m = 
Pto = 1 and write: 

A{w-w') = j djdf Py,,n,'{j,j')c{a^{j))c{ay,'{j')) (77) 



1. Calculation ofA{w) for class I 

We define ai := ayj{j), 02 := a^'(j), 03 := ««;'(/)• 
Given the previous remark to compute the cumulants we can 
set 

j — w = lnai , j — w' = ln 02 (78) 
j'-w'=\na3 , 02/01 =e~^, (79) 

where we denote W = w' — w > 0. The joint probabiUty (71) 
then reads 

Pw;w'{j,j')djdj' 

= dai{l - c-'^)da3e-"i(i-<'""')-''36»(/ > j) 

- j)dj'e-^daie-''^ . (80) 

This yields the second moment 

{w-j){w'-j'f = [1 - e-^] / dai / das 

Jo Jaic-W 

X exp (— oi [1 — e~^] — as) In ai In 03 



+e ^ {Ww — j + (w — — w — j 



(81) 



Note the integration interval for 03 which corresponds to j' > 
j using (79). We recall that 



w — J = — I dae \na = 
'0 



(82) 



-3?= / 
Jo 



da e-"(ln af =-(%+'— . (83) 



Thus we obtain: 



A{W) = (1 - e-^) / daie-"i(i-^""') Inai 
^0 



roc 
J ai i 



das e In 03 



+6-^^(71 + -^ + W^7s)-7l (84) 

The calculation is performed in Appendix B. The final result 
for the fixed-point function of class I is: 

2 2 
AH = ^+Li2(l-e'") + ^ (85) 
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where Li„(z) = Yl'kLi -^'^/fc"- One can also use the alterna- 
tive formula (B5). Another equivalent compact expression for 
the result is: 



AH 



El + nw 
5 — ( 

n=l 



(86) 



A numerical test is performed in section V, see figure 1 1 . 
The behavior of the fixed-point function at small «; > is: 



2 ^ ^ 7 

r , , 77" W W W W 

Aiw) = w-\ 1 

^ ' 6 4 36 3600 211680 



hence we confirm that there is a cusp, and a power series ex- 
pansion in I . The behavior at large w is easier obtained from 
(86) and reads: 

K{w) = {w + l)e-^ + ^(2w + l)e-2«' + \{;iw + Ije-^™ 
+0 (e-^"") (87) 

It is characteristic of short-ranged correlations in the force 
with an exponential decay. 



2. Calculation ofA{w) for class III 

As for class I, we define ai = aw{j), 02 = a-w'ij), 0,3 
dw' W) with (attention: j is ahead of w) 



j -w 


1/a 


(88) 


j - w' 


1/a 
= 0"! 


(89) 


j'-w' 


1/a 
= % • 


(90) 



One must distinguish the cases j > w' (then f > j is equiva- 
lent to as > a2 and W = w' — w = o}^" — Og^") and j < w', 
in which case a2 = 0. There are thus two pieces for the part 
with / > j: 



di := (j - w){j' - w%'>j 

fw' J l-CO 

= / di^e-«i / da3e-«3(aia3)* 
Jw Jo 

Jw' dj 



(91) 



In the first integral the integration bounds over ai is from 



to W" and in the second the relation w' — w = aj^" — al/" 



holds. Hence one obtains 

^ r-W 



di = r ( 1 -h - 1 / dai a^^e-"-' 



a) Jo 

7v^o V dai J 
The second contribution {j = j') is: 



(92) 



a 



d2 := (j - w){j' - w%'=j 
Jw dai 



(93) 



One has for the disorder correlator: 



A(W) =di+d2-[ I daie-"-'a/' 



(94) 



Using 02 = {a/" - W)°', daz/dai = (1 - wa'^^"')"'-'^ and 
the variable y = a^^", a^" = y + W, one finds: 



A(W) = -T (1 + i )r( l + 

\ a J \ a 

/•OO 

+a / 6.y{y + W)e-^y+'^'^°' 
Jo 



(95) 



X [y« + ey"T{l + -,yn {{y + W)"-' - y"-^) ] 

a 

Integration by part of the last term finally yields: 

AH = -r(i + i)r(i + l,.^ 

+wT[ l + i'je-"'" 



+ [ d2/e-(''+"')"+2'°r (1 + 
Jo \ a 



We recall that 



\a,x)= f 

J X 



(96) 



(97) 



Hence we find a fixed-point function continuously dependent 
on a, and C, = 2 — 2/a yielding a unique form for each value 
ofO < C < 2. 
The value at tz; = has a simple expression: 



A(o) = r(i + -) - r(i + l 

a J \ a 



We find that the function A{w) has a cusp with 



-A'(0+) = 



a 



(98) 



(99) 



Since there are a priori w°' terms (we recall a > 1), we want to 
understand at which order the expansion in \w\ breaks down. 

More explicit expressions can be obtained in special cases. 
For the box-distribution, i.e. a = 2, we find 



A{w) = 2w-e'" y/iT {2v? + l) erfc(w;) + 



1 r- 
+ 2^ 



w e 



(100) 



which has a power series expansion in around w = 0: 

/ttu;^ w'^ y/nw'^ 
'24 30"'' 120 



r , . I ^ TT\ WTTW W 

^(") = (i-4)-V + T- 



'^210 672 1512 ' 4320 ' 11880 
+0(u;") (101) 



11 



(100) together with (70) is checked numerically in section V, 
see figure 1 1 . 

However for a non-integer the situation is more compli- 
cated. Despite the presence of terms (we recall a > 1), 
the one sided second derivative at A"(0+) seems to exists. 
Numerically one finds that for a non-integer, but close to 1, 
e.g. a = 3/2, the third derivative at exists, but does not go 
to zero with a finite slope. Analytically, one obtains an expan- 
sion around a = 1, as 

A{w) = Ao{w) + {a- l)Ai(w) + 0{a - if (102) 

with 

AoH=e-"' (103) 

Ai{w) = {w - 2)T{0,w) + e""' (- {w + 2) log(w) - 4) 

(104) 

with the incomplete F-function defined in (97). Note that 

dl AAw) = -e-'-iwlogiw) + 1) 
and the Taylor expansion of (104) around is then: 



Ai(w) = (-4 + 27) + {2- 7)w - 



w 



+ 



11 

36 ~ 

+ O (w^) 



log(w) 



log(w) 
12 



13 
144 



whose logarithmic part can be summed up as: 

Ai(«;) = [2-w-e-'"{w + 2)] log(M;) 
+ analytic in |w| 



(105) 



(106) 



These expansions confirms that for the Weibull class at non- 
integer a the second derivative at 0+, A"(0+), exists, but not 
the third one. 

Finally let us note that if pj^^ is scaled as a one 
should recover class 1 from the large-a limit of class 
III, i.e. one can indeed check from (79) and (109) that 
lima^oo a^A''''"(w/a) = A'(w). (the indices refer to the 
class). An example of these limit procedure is given below. 



3. Calculation of A{w)for class II (Frechet) 

We define again ai = a^ij), 012 = a^'ij), as = a^'if) 
with 



j -w = 


-l/a 


(107) 


j -w' = 


-l/a 


(108) 


f = 


-l/a 
— 03 


(109) 



One sees that j can only vary in the interval [—00,10], and that 
j' > j is equivalent to as > 02- As for Weibull, there are two 
pieces: 



di := {j - w){j' - w')y>j 



(110) 



, d(ai-a2) ^_(„,_„,) rdase-«3(aia3)^ 
J 02 



where the relation W = w' — w = a2^^°' — a^ holds. This 
yields: 



di= I dai{l-^] aj"'/"e-«i+»^r(l - i as) (111) 



da2 

/o \ dai 

The second contribution is: 



d2 ■■= (j -w)ij' -W%>^j 

-I 

-L 



di^e--(aia2)-i/° 
da2 



dai-^e-"Uaia2)-^/" (112) 
dai ' ' 



One has for the disorder correlator: 



A{W) =di+d2- i^j daie-"ia^^/"J (113) 
Using the variable y = a^^^", a2^^" = y + W, one finds: 

POO 

A{W) =a dyy{y + Wye-^" 
Jo 

POO 

+a dyyiy-'-'^-iy + W)-'-") 
Jo 

-y-'+{y+wr''T{l -^,{y + W)-") 

(114) 



X e 

1 \2 



As for the Weibull class, see (96), we can integrate by part the 
second term into: 

A 

-^lim^ J dye-^""+(^+^)"°^ [yr(l - i (y + H^)"")] 


+ Ar(i-i), (115) 

where the last term comes from the upper bound in the partial 
integration, after putting y ^ A ^ 00 there. Rewriting it 
as dj/ r(l — l/a), and using the fact that the first term in 
(1 14) cancels, we arrive at the simple final expression for the 
Frechet class: 

AH = -r(i-i)2 (116) 

00 

+ 1 dy r(i - 1) - e-^""+(f+-)""r(i - i (2/ + w)-") 



We find a fixed-point function, which depends continuously 
on a, hence on ( = 2 + 2/a with a unique form for each 
value of C, in the domain < C < 2. 

We can obtain the small-w behavior easily from (1 14). The 
value at zero. 



A(o) = ^(l-^)-^(l-i)^ 



(117) 



is found consistent with the results from the previous section 
on the distribution of the critical force. 



(118) 
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One sees that A(0) is defined for a > 2 and diverges as a ^ 
2+. As discussed in the previous section this is because in the 
Frechet class the distribution of the critical force has algebraic 
tails and an infinite n-th moment for a < n. We thus consider 
a > 2. The fixed-point function has a cusp, with from (116) 



_A'(0+) = ^ , 

a 

and a well defined Taylor expansion in \w\: 

u2 rr-- 



(119) 



AH^r(H-2)_r(i-i)2 + 

(a + l)r(l)i.3 



aw 
Aa + 



36a + 24 



48(2a(4a + 9) + 9) 



a^(a + l)(«^-4)r(2+|) 
240(15^2 + 32a +16) 



+0 (w^ 



(120) 



The large-w behavior of the fixed-point function is quite dif- 
ferent from the other classes. Indeed A(w) decays to zero 
rather slowly for large w. This can be seen by writing 



fiy+w) 

ti{w)= I dye"^"°+(^+"')"" 



1 



(121) 







The leading term for large w comes from ti{w), via a series 
of approximations: 



dye 



-y °+(y+«') 



(y+w) 



dtr 



a-l 



a- 1 



dy {w + yy-" 



(a-2)(a-l) ' 



(122) 



whereas t2{'w) is of order 1/w" plus a constant, since one can 
simply expand the exponential function for large w. For a — 
3 the function A{w) and the asymptotics (122) are plotted on 
figure 6. 



IV. AVALANCHE-SIZE AND WAITING-TIME 
DISTRIBUTIONS 

A. Avalanche-size distribution 

Successive avalanches, or jumps, occur at a discrete set of 
w = Wi such that u{w~) = ji and u{w^) = ji+Si where s = 
Si > is the size of the avalanche. The waiting time between 




FIG. 6: A(to) from eq. (116) for a = 3 (bold). One clearly sees the 
long tail for w ^ oo. The asymptotic behavior for large w from eq. 
(122) is shown (dashed/red) as the small-ui expansion (blue/dashed). 



F(u}k 




FIG. 7: Geometrical construction for the size of an avalanche 



consecutive avalanches is denoted W = Wi — Wi-i. It is not 
properly a time, but we will term it here loosely waiting time 
since for a driving with a constant velocity it is the waiting 
time = W/v, with here v ~ (in that limit the jump time 
is negligibly shorter). Here we compute the joint distribution 
of avalanche sizes and waiting times. In section V we discuss 
an algorithm to generate the sequence of avalanches and the 
Markov-chain property. 

There are two useful probabilities for which the general ex- 
pressions for the discrete model are easy to write. The first is, 
for w' > w. 

P^{r,w'):= (123) 
Proba(u(?«) — j and next avalanche is in [w' , oo]) 

When u{w) — j the next avalanche occurs at w" > w such 
that m^{j — w") = Fj. Thus to realize w" > w' we need 



Fj < m {j — w'). Hence one has: 



PUj;w')^{l-Him\j-w')) II Him^ik-w)). 

(124) 

For w = w', Pw{i\ w) = Pw{j) defined in (17). Of interest is 
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the differential waiting time distribution 

oo 

PUw')dw' = -dw'd^, J2 Pw{j;w') (125) 

j=-oo 

= Proba(given w, the next avalanche is in [w', w' + 6.w']) 

The joint (integrated) probability is more involved: 

P^{j,w',S)dw' := 

Proba(u(w) = j and next avalanche is in [w' ,w' + dw'] 
and of size s > S ) (126) 

One has: 

s 

P^{j,w\S) = n H{m\p + j - w')){-d^,P^{j-w')) 
p=i 

(127) 

and of course PwiJ^w^S = 0) = —dw'Pw{j',w'), i.e. the 
probability that u{w) = j and that the next avalanche occurs 
at w' (per unit dw'). Of particular interest is: 



PUw',s>S) :=^P„0>',5) 



(128) 



= Proba(next avalanche is in [w' , w' + dw'] 
and of size s > S) 

This is illustrated on figure 7. 

We now consider the limit m ^ in which the contin- 
uum limit can be used, and extract the distributions of waiting 
times and avalanche sizes. We find that they depend on the 
same single scale pm as defined in the previous Sections. The 
waiting time distribution is denoted below P{w) (and some- 
times P{W)), and has thus the scaling form: 



Piw) = p-'P{w/pm) 



(129) 



However, to simplify notations, unless specified, we drop be- 
low the tilde on P and formally replace pm 1 (formally, 
since all statements below are about scaling forms in the linut 
m ^ 0). When no confusion is possible we use the sym- 
bol w for either the argument of u{w) and the waiting time 
w ^ Wi — u>i_i, and use W for the waiting time when con- 
fusion is possible. Similarly the avalanche size distribution is 
loosely noted with the same symbol: 



(130) 



and we perform the same simplification in notations, and sim- 
ilarly for the joint distribution. The dependence on pm can be 
restored by replacing in the final formulas 



J 



w 

Pm 



(131) 

(132) 



similarly for S and W, and correcting the normaUzations of 
probabilities to 1. 
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FIG. 8: The probability Pw{w') to have a point w in the interval of 
size W preceding w'. 



B. Waiting-time distribution 

Using the method introduced in the previous section, see 
especially Eq. (21) ff., one finds that the continuum hmit for 
the probabiUty Pw{j', w') is 

P^(j- w')dj = da„/(,7)e-''™(^') , (133) 

and we use again the notation ai = a^ij) and = aw'{j)- 

1. Class I ( Gumbel) 
We recall (75) for class I 



j -w = -/c -l-lnai 
j-w' = -f^+\na2 



(134) 



where /° is a non-fluctuating constant, and as in (79) a2 
e~^ai. Thus one finds: 



P^{j:w'}dj = e-^daie-''' 
Integrating over ai in [0, oo[ one obtains 

P^{w') = e{w' -w)e'^-'^' . 



(135) 



(136) 



From (136) we can infer the distribution P{W) of wait- 
ing times W. Since the probability that a uniformly cho- 
sen w on the real axis falls in an interval of size W is 
WP{W)/ JdWWP{W) (see figure 8), and that then the 
probabiUty of of w' — w is uniform, i.e. 9{0 < w' — w < 
W)/W, multiplying the two expressions one finds the gen- 
eral relation: 



n dWPiW) 
dWWPi^) 



(137) 



This is true only in the small-m hmit (because of uniform 
measure assumption). For any member of class I we find the 
distribution of "waiting time" W (restoring the dependence 

on Pm): 



(138) 



2. Class III (Weibull) 



Let us consider now class III. One has, foUowing the nota- 
tion in section III D 2: 



J - w = -f^ + a/ 



J - w 



(139) 
(140) 
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Using the variable y > such that a^^" = y and a^" = 
y + W one now finds from (133): 

PUj;w')dj = ay"-ie-fe+>^)° dy (141) 

Integrating w.r.t. y, taking —dw' and integrating by part we 
obtain: 

/>oo 

P^{w')=a{a-\) d2/y"-2e-(^+"''-"')" (142) 

We recall that a > 1 hence the integral is always convergent. 
From (137) we obtain: 



{W)-^ =P^{w)=aT[2 



1 



a 



(143) 



and, taking —dyji on both sides of (137) we obtain for class 
III: 



P{W) 



r dyy'^-^y + Wr-'e-^y+^r 



(144) 



If a > 2 it can be integrated by part into 

P{W)=^^^§^ r dyy-^e-^y-^r _ (145) 



For completeness let us give the result for a = 2, which cor- 
responds to a box distribution for the forces: 



P^{w') = V^erfc(iy) 
which is a simple one-sided Gaussian. 



(146) 
(147) 



(148) 
(149) 



3. Class II ( Frechet) 
Let us consider now class II (Frechet). One has: 

-l/a 

J - w = -a^ 

I -l/a 
J -W = 

Using the variable y > such that a^^^" = y and 03 = 
y + W one now finds from (133): 

PwU; w')dj = a{y + W)-"-^e-y~''dy (150) 
Integrating w.r.t. y, taking —8^' we obtain: 



Py,{w') =a{a + l) 



iy + w) 



2+a 



(151) 



We recall that a > hence the integral is always convergent. 
From (137) we get 

(W^)-i = P^{w) = ar(^2 + . (152) 

Deriving again (137) w.r.t. W, we obtain 



P{W) 



(a + l)(a + 2) 



dy 



r(2 + i) io {y + wf+' 

One easily checks that (153) is correctly normalized. 



e-^' . (153) 



C. Joint avalanche-size and waiting-time distribution 

The continuum Umit of Eq. (127) can be written as 

P^{j, w', S) dj = [-d^>{da2 e-'^i)]e-"^+»= , (154) 

where we use the same notations (j) = ai, ciw' (i) = 0-2 as 
above and in addition a„' (j' = j + S) — as. We remind that 
in order to include Pm, S gets rescaled as j and w. 

I. Class I ( Guinbel } 

Consider now class I. In addition to (134) one has 

j + S-w' = -f° + lna3. (155) 
Hence we have 02 = e~^a\ and 03 = e~^~^ai where 

P^{3,w',S)dj = [-a^,e-^(daie-''0]e-"^<^""('''-^) 

= e-^dai exp ( - ai[l + e-^(e^ - 1)]) 

(156) 

Integrating over j and ai respectively, we obtain the joint dis- 
tribution: 



P^{w',s> S) = 9{w' -w) 



1 



+ e^' - 1 



(157) 



Again this is the probabihty that if one observes the system 
at w, the next avalanche occurs at w' and has size s > S. 
This allows to find the joint probability P+ [W, s > S) that 
Wn+i — Wn = W and the next avalanche is s„+i > S (see 
Appendix V for further notations and definitions). Indeed one 
has, by the same reasoning as above [94]: 



C_^dWP+{W,s>S) 
dWWP+(W,0) 
This yields for class I: 

P+{W,s> S) = 



= Py,{w' , s > S) (158) 



{e^ + - 1)2 



(159) 



Integrating over W this yields the cumulative joint waiting- 
time and avalanche- size distribution for class I: 



P{w >W,s>S) = 



1 



exp(VK) + cxp(5) - 1 ■ 



(160) 



Setting W = gives the (rescaled) avalanche-size distribution 
for class I: 



P{s) = e-' , P{s >S) = e- 



(161) 



Setting S = 0, and deriving w.r.t. W reproduces the waiting- 
time distribution (138). We thus obtain that the avalanche ex- 
ponent, such that P(s) ~ s~'^ at small 1 ^ s <C p™, is here 
r = 0. We can give the lowest moments: 



(ws) 



1.64493 



(w) = {s) = l 



(162) 

(163) 
(164) 
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in units of pm- We also note that the relation: 



is) 



-2A'(0+) 



(165) 



is obeyed, using A'(0+) = — 1 from (87). This relation be- 
tween the cusp and the second moment holds quite generally 
[62J and is used here as a useful check. 



2. Class III (Weibull) 
Consider now class III, i.e. (140) and 



(166) 



This leads to 

Pw{j,w',S)dj 



xe 



(167) 



Now setting as before y = Oj^" = j — w' , the only derivative 
in the bracket to be taken is of the variable y, i.e. —dyj/y = 1. 
In short-hand then gives 

P^{j,w',S)dj = -a^,ay"-Mj/e-(^+^)"-(f+^)"+^" 

(168) 

Integrating over y (i.e. j) it yields 



/■OO 

P^{w',s> S) = a{a-1) / dyy' 
Jo 



,a-2g-(j/+W)°-(3/+S)°+j/° 



(169) 

The final result for the joint probability for class 111 takes the 
form: 



P+{W,s> S) 

pCX) 



(170) 

a-2- {y+W)''-{y+Sr+y'^ 



Integrating over W yields the (rescaled) cumulative distribu- 
tion for class 111: 



P{w > W, s > ,S) 
a 



(171) 



a-2 -{y+W)'='-{y+S)"+y° 



dyy'^-'e 



r(i - i) Jo 

Setting = 0, we obtain the avalanche- size distribution. 



P{s > S) 



r(i - i) Jo 



,c-2-{y+Sr 



dyy"-'e 



(172) 



On the other hand, setting S = in (171) gives back the 
waiting-time distribution (144). 

Let us comment these results. Since (172) can be Taylor 
expanded in S around 5* = 0, it is clear that for class 111 also 



the avalanche exponent is again r = 0. For large S, the decay 
of (172) is P{s > S) ^ P{S) ~ exp(-5"), i.e. a stretched 
exponent decay with exponent 6 — a. Next, one checks that 
the general relation involving the cusp is obeyed: 



(^^)_ r(2 + i) 

(s) a 



2A'(0+) . (173) 



using (99), a useful check on our calculations. 

One can also obtain simple expressions for the lowest mo- 
ments (in units of p^); 



is) = H = 

{sw) 



1 



"r(2-i) 
r(i + i)if_i 



r(-^) 



(174) 
(175) 



as a function of the harmonic number Hn = J2k=i 1/^' 
H_L :=7E-hV(l-l/a),V'(a;) := T' (a;) /F (a;). From these 
one can construct a fuUy universal dimensionless ratio: 



{sw) 



"'sin (I)" 



(176) 



Finally one can check that for a — > oo, one recovers class I 
distribution (160). More precisely: 



lim P{w > W/a, s > S/a) 

a. — *oo 

1 



(177) 



= lim 



oor(2-i) 



dx X 



X exp —X 



1 + 



W 



s 



- 1 



/•OO 

= dx exp {-X [e^ + e"^ - l] ) = (160) , (178) 

where we have used x = y" as variable and lima_>oo(l + 
W/a)" = exp(W^) and that the factors of x^^" could be 
dropped. 



3. Class II ( Frechet) 
Consider now class 11, i.e. (140) and 



-l/a 



(179) 



j + S -w' = -Pmds 

hence S = a2^^°' — a^^^" = y + W — a^'^'"^. This leads to: 



-l/a 



P^{j,w',S)dj = [-d^,aiy + WT''-'dye-y~^] 



X e 



{y+W)-'='-iy+W-S)- 



= a{a + l){y + W)-''-^dy e'^ 



X e 



{y+wy-iy+w-sy 



(180) 
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The domain of integration depends on whether W > S or 
S > W, and can be expressed as 



> S) 
= a{a + 1) 

— a{a + 1) 



dy 



X e 

oo 



-y-'' + {y+W)-''-{y+W-S)-'' (•^g^-j 

dy 



m&K{W,S) y 



Q + 2 



(182) 



This translates with the same arguments as for (171) into 

P{w >W,s> S) 

OO 



(183) 



dy 



.a+2 




FIG. 9: The plot shows u{w) for the forward motion with 
(w — u) = fc, and the backward motion with (u — w) = fc. 
Both trajectories are indicated with arrows. The world line u{w) 
fluctuates around w ± f^m^'^ . These fluctuations have a geometri- 
cal interpretation as the yellow/shaded region, whose area is related 
to the expectation (ws). 



As for the Weibull and Gumbel universality classes, this ex- 
pression is symmetric in S and W . We therefore conclude that 
for all three classes: 



P{w) = P{s) 



(184) 



This property is proved with slightly more general argument 
in Appendix D. It is of course valid only for the rescaled 
distributions in the limit to ^ in the sense described above 
(at the level of the discrete model is a continuous variable 
while s is discrete). 

The resulting distribution of (rescaled) avalanche size thus 
reads: 



P{s) 



r(2 + i) 

with (in units of pm)'- 



dy 



{y + s) 



= ar( 2 + - 

a 



-2r(-l/a) 



= -2A'(0+) 



(185) 

(186) 
(187) 



the last equality being a check, using (119). Note that P{s) 
has power law decay for large s, i.e. P{s) w 5^^^+"-'. The 
avalanche exponent however is still t = 0, since it is related 
to small avalanches (i.e. s <C Pm). The Frechet class yields to 
a cutoff for large avalanches (i.e. s ^ p„i) which is itself a 
power law. 



a symmetric distribution P{~f) = P{f)- The plot shows 
u{w) for the forward motion with m? {w ~ u) ~ fc > 
(disorder or translational averages) and the backward mo- 
tion {u ~ w) — fc- Both curves u{w) fluctuate around 
w ± fci'm?, and the enclosed area A of the hysteresis loop 
exhibits a uniform part A = 2|/"(to)|w, computed in Section 
III C for each class, plus a fluctuating part A: 



A^A + A 



A 



WiSi 



(188) 



It has a geometrical interpretation as the area of the yel- 
low/shaded region. Its average value per unit length over a 
large sample w e [0, AI] is 



1 f\) - (WS) 



(189) 



since the number of avalanches is — J2i — 
hence (s) = M/N. Note en passant that one has also 
~ M hence (s) = (w), i.e. the first moments of s 
and w are always equal. (We neglect boundary contributions 
which for uncorrected disorder scale subdominantly). Sim- 
ilarly one can consider the moments of the local hysteresis 
area: 



(190) 



4. Local fluctuations of the area of the hysteresis loop 

Figure 9 illustrates a typical hysteresis loop for uncorrected 
random force landscape, for convenience assumed to have 



They can be obtained from the moments of the variable a = 
ws. Hence it is useful to compute the distribution of this vari- 
able for the three classes. This is performed in Appendix D. 
Let us give the result for the (rescaled) distribution of the 
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A(«) 
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FIG. 10: Numerically calculated A(u) as a function of u (red, fat) 
for the box distribution P{f) = Q{f)0{l — /), and comparison 
with analytical result (black, dashed) from equation (100); A — 1, 
a = 0, A = 1/2, pm = \/2/m, = 10^^; there is no adjustable 
parameter. 



Gumbel class: 



oo 



P{a > A) = 




(191) 



from which we give some moments: 

(a) = 1.64493 , (a^) = 18.3995 , 

(a^) = 547.343 , {a'^) = 30764.6 , (192) 

measured in units of . 

V. NUMERICS FOR THE TOY MODEL 

A. Basic definitions 

In this section we study numerically the discrete model of 
the last section. We consider four different disorder distribu- 
tions: 
Box : 

Pi(/) = lif/e [0,1], Oandelse. (193) 

Exponential: 

P2 (/) = cxp(/) if /> and else. (194) 

RB (random bond): 

short-ranged correlated potential. The resulting force at site 
i is, fi = Ci — Ci+i where the energies are uncorrected 




FIG. 11: A(m) as a function of u for the exponential distribution 
(red, fat), and comparison with analytical result from equation (85) 
(black, dashed); — = 0.003, no adjustable parameter. 
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FIG. 12: A(ii) versus « for a = 3 (Weibull class); = 0.0001. 
No adjustable parameter. Red = data with error-bars. Blue = an- 
alytic solution (96) for a = 3. The other parameters used in (50) 
ff. are: A — I/7 and a = ^ — 1. We have chosen 7=1/2 
as an example, hence a = 1, a = A = 2, A = ^. We use 
= 0.0001, with 10^ disorder points. One finds /c/m^: -599.8 
(numerics) versus —597.79 (analytic). The parameter-free numeri- 
cal result for A(u) is compared with the analytical one, and found in 
excellent agreement. 

random variables distributed with the box distribution Pi(e). 
We call this distribution Psif). 
Class III with a — 3 : 

P4(/)-2/e(0</<l). (195) 

Note that power law-distributed forces can be generated by 
defining / :— x'', with x £ [0,1] uniformly distributed. This 
yields 

^^a) = -/-~'0(O < / < 1) , (196) 

7 

with 7 = 1/2 for (195). 

We integrate numerically the equation of movement (8), 
first at 77 = 0. This is the discrete model defined in section 
III A. In practice, for given w, we move the particle as long as 
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P{W),PiS) 

1.2 - 



Pif) = e(o < / < 1) 




w,s 



FIG. 13: Reseated disorder correlator: The function Y{u) is the 
rescaled version of A (it), s.t. Y{0) = 1 and Y{u) = 1. The 
data is the same as on figures 10, 11 with in addition RB disorder 
(i.e. random potential): Blue is the box distribution in [0,1], red the 
exponential, green is RB. This shows that for m — 0.003 the differ- 
ent microscopic disorders yield very similar rescaled correlators, al- 
though the unrescaled ones are different. For the difference between 
rescaled disorders see figure 14. 




FIG. 14: Comparison of the analytical results for the various classes. 
The function Y{u) is the rescaled version of A(ii), s.t. Y{0) = 1 
and J^^Q Y{u) — 1. Plotted here are the (small) differences ^1(2;) — 
^111,2(2;) (red) and 1^1(3;) — 1^11,3(2:) (blue). The index refers to class, 
and a. This explains why the various classes are very close on figure 
13. 



the force, i.e. the r.h.s. of equation (8) is positive. The point at 
which we stop defines u{w). We then update w ^ w + 1. This 
is an approximation to the process defined in section III A, but 
since jumps as well as waiting times diverge when m 0, 
the scaling limit is the same. A second algorithm, described 
in Appendix C, was used to independently compute P{w) and 
P{s) (not shown), and check the present results. 



FIG. 15: P{W) (orange dots) and P{S) (green dots) as a func- 
tion of W and S respectively, for the box distribution P{f) = 
0(/)6(l — /), and comparison with analytical result (black, dashed) 
from equation (100); A = 1, a = 0, A = 1/2, = \/2/m, 
— 10^^; there is no adjustable parameter. One sees that the two 
distributions are identical, with P{W) plotted on top of P{S). How- 
ever for small avalanche sizes S, the numerics has not yet converged 
against the analytical result, i.e. finite-size corrections are present and 
visible on the plot. 



adjustable parameter in figures 10 to 12 and the agreement is 
excellent. 

An general important question at depinning is whether the 
random bond class (i.e. uncorrected potentials) flows to the 
random field one (uncorrected forces). This appears clearly 
in figure 13, where we plot the rescaled (as explained in the 
caption) A(u), for the 3 disorders P2{f), and Pzif), 

defined in (193) ff. The cross-over as the mass decreases from 
RB disorder to RF disorder can also nicely be seen in our sim- 
ulations, presented on figure 2. One expects that the random 
short-ranged energy model, i.e. fi — Ci — e^+i with dis- 
tributed with Prb (e) should flow to the random force model 
with P{f) — /g -Prb(/ + e)PaB(e), i.e. the convolution of 
-FRB(e) and Prb(— e). This is because the rare large forces 
are isolated and become uncorrected. This predicts that the 
box distribution for e should flow to the a ~ 3 class III. 

We have seen on Fig. 13 that the shapes of the correlator 
functions A(u) (i.e. their rescaled form as explained in the 
caption) are rather similar for the various universality classes, 
while their unrescaled forms are very different. These rescaled 
forms obtained from the analytical calculations are compared 
in Fig. 14. 



C. Shocks and Avalanches 



B. A(it;) 

We have shown on figures 10, 11, and 12, comparisons be- 
tween the numerically computed functions A (it), and the an- 
alytical predictions. The corresponding analytical results are 
referenced in the corresponding captions. Hence there is no 



In figure 1 we have shown the avalanches, also called dy- 
namical shocks: As a function of w, we plot w — Uw (minus 
its average), for different masses. First consider the smallest 
mass, — 0.03. We see that w — is growing linearly 
with w, before it jumps. The linear parts are those, at which 
the particle is localized by a large disorder force, before jump- 
ing (in zero time at rjv = 0) to a new position (vertical parts). 
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When decreasing the mass, we see that the linear parts, i.e. 
the "time" (i.e. distance in w) between jumps, as well as the 
jumps itself become larger, while sharing parts of their trajec- 
tories. This can be interpreted as merging of the (dynamical) 
shocks. 

In figure 15 we show the avalanche-size and waiting-time 
distributions for P{f) being the box distribution. It clearly 
shows that avalanche-size and waiting-times distributions are 
identical, as follows for all disorder classes from (D3). More- 
over, the result is in agreement with the parameter free predic- 
tion of eq. (172), using pm defined in equation (59). 

D. Finite velocity 

In this section, we consider the equation of motion at finite 
velocity v, or rather at finite rjv, since the latter is the parame- 
ter entering all equations. 

The algorithm for finite rjv is as follows: We generalize 
the position ut = u{wt) of the point to now take non-integer 
values. It follows the Langevin equation 

V^ut = /[„,] + m^{ut - wt) (197) 

where 

[u] := largest integer < u . (198) 

Condition (198) reflects that the disorder only changes at dis- 
crete values of u. 

In practice, we discretize( 197) with a step-size (5f — 1/100, 
integrating this discretized equation of motion in time follow- 
ing the Ito scheme, and using wt = vt. 

V [ut+st - Ut] = St + m^(ut - Wt)] (199) 

To guarantee that our time discretization is fine enough, we 
report "maxslide", the maximum of (199) encountered in a 
simulation. For the simulation shown on figure 16, this was 
0.05. The figure shows A(u) at j]v — 0, calculated with the 
discrete algorithm used in the previous sections, and A(u) at 
?/ = 1, and V = 0.2. The microscopic disorder is a box dis- 
tribution for the force, given by eq. (193). The result clearly 
shows a rounding of the cusp by the non-zero velocity. 

VI. LONG-RANGE CORRELATED FORCES 

In addition to the three universality classes for short-range 
(SR) correlated forces, there is also a family of univer- 
sality classes for long-range (LR) correlated forces. Con- 
sider a gaussian distributed force landscape with no bias 
F{u) — F{u') — and second moment: 

[F{u) - F{u')Y = 2a\u~u'y (200) 

We focus on 7 = 1, i.e. a Brownian-force landscape, but one 
expects a continuously varying fixed point as a function of 7. 
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FIG. 16: Rounding of A(it) through a finite velocity; rjv = 0.2. 
Inset: A'(ti). The disorder distribution is P(/) = 9(0 < / < 1), 
= 0.003. 



Although in most cases random-force landscapes at depinning 
have short-ranged correlations, these more exotic LR land- 
scapes exhibit some interesting properties. Note that in the 
statics the (very) LR correlated random potential landscape 
corresponding to the case 7 = 1 was studied by Sinai [71]. It 
was found that shocks are dense there. 

Quite remarkably, exact results can be obtained for this 
model for any non-zero velocity v > 0, as first noticed by 
Alessandro, Beatrice, Bertotti and Montorsi [77] who intro- 
duced this model, hence referred to as ABBM model, as a 
realistic description of the Barkhausen effect in metallic ferro- 
magnets, and compared the results to experiments. The mass 
term originates from the magnetostatic fields: the demagne- 
tizing field (resulting from the effective monopoles sitting at 
the end of the sample) provides a long-ranged restoring force 
which acts as a spring, precisely as in the model considered in 
this paper (see also discussions in [79] and [78]). 

We first obtain some results on the quasi-static model, then 
we recall some results of the ABBM analysis at v > 0, and 
obtain from it the renormalized correlator at non-zero velocity. 




FIG. 17: defined in (201) for /a random walk with D = 1.67x 
10^ (blue). The dashed (red) fine denotes ^('u) = u. We show 
expficitly the graphical construction of u(w) for w = 3400, as well 
as a larger avalanche of size S = 3785. The system is discretized 
with stepsize 1 in ti direction. 
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We then discuss how these results match at t; — > 0+. 



2. 2-point conditional distribution 



A. Quasi-static motion 



The forward process u 
root of 



is defined as the smallest 



(j){u) = u — m ^F{u) = w . 



(201) 



For 7 = 1 the process (p{u} is a Brownian motion (BM) of 
diffusion constant D = a/rn^ and upward drift b = 1. The 
time of the Brownian ist = u. From Fig. 17 one sees that it 
becomes a first-passage-time problem, i.e. u = u{w) is the 
first time the process </>(■«) reaches altitude w. Conversely 
the process w{u) = max„/<„ (f){u'), is the maximum posi- 
tion over all previous times reached by a brownian motion. 
To avoid pathologies we assume a cutoff which makes F{u) 
smooth at very short scales; equivalently we can discretize in 
u-direction, as was done to generate figure 17. In Appendix 
E we have collected some useful properties of first-passage 
times and maxima of the Brownian motion that we now use 
extensively. We refer to this Appendix for all details. 



1. Avalanche distribution 

From Fig. 17 one sees that the avalanche distribution P{s) 
identifies with the return probability to the origin of the BM 
with a drift h=l and diffusion constant D = <j/m^. 



P{s) « Wo) 



1/2 



s ^/^exp 
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(202) 



where u — s, P{u; W) is defined in (E7), 6=1 and D = 
a/m'^. Wo is a non-universal short-distance scale. We have 
defined sq ■= Wq /D and Sm '■= D the short-scale and large- 
scale cutoffs for the avalanche size. In the limit of small m 
one has sq "C Sm. allowing to drop the last term in (202). 
There are many small avalanches of the order of sq, i.e. the 
distribution is concentrated at sq. However the moments (s'p) 
for p > 1/2 are dominated by large avalanches. For s ^ sq 
one has: 
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(204) 



and (s) = y^SmSo- This is exactly the distribution found in the 
mean-field theory of sandpiles [80] and of the random-field 
Ising model, and was recently shown to hold, using FRG [62, 
84], for elastic manifolds in d = 4. The random-walk picture 
goes back to the so-called Gallon process [83] for survival of 
family names (see [81] for a recent discussion in the context 
of depinning) which exhibits the same mean-field power-law 
behavior at the threshold. 



It turns out that the 1-point-probabihty and critical-force 
distribution is a subtle issue for this model, due to the long- 
range nature of the landscape and the choice of boundary con- 
ditions. We do not discuss it in details here, but some consid- 
erations are given in the Appendix F. A full solution requires 
a separate study. 

We can still offer some simple remarks. If we know, e.g. 
by observation in a numerical simulation or an experiment, 
that the process is such that u{wi) = wi, then one can easily 
compute, from the Markov property of Brownian motion, the 
probabilities of a\\ future events, i.e. the conditional probabil- 
ity for u{w2), ■ • ■ , u{wn) with W\ < W2 < ■ ■ ■ < Wn'. 



XU2, ..,Un\u{wi) = Ui) 



(205) 



P{U2 -Ui;W2-Wi)x ...X P{Un - «„_!; W„ - Wn-l) 

where P{u; w) is the first-passage-time probability defined in 
(E7). Computing the moments one finds: 



u{Wn) -Wn= U{wi) - Wi 



[u{Wn) - Wn - u{Wp) - Wp]'^ = 2D\Wn - Wp\ 



2a 

= :^\^n-wp 



(206) 



(207) 



for any n,p> 1. This defines the renormaUzed correlator: 

A(0) - A(w) = aw , (208) 

i.e. it is exactly the bare disorder correlator. Note that al- 
though this result was derived from a conditional probabil- 
ity, it is independent of the choice of ui and wi provided all 
points Wi to which it applies are larger than wi . Note also 
that although the 2-point correlator is the bare one, the higher 
cumulants are different: they are non-trivial and can be easily 
computed from (206). Hence the system flows to a non- trivial 
fixed point. Remarkably, one can check that (208) is an ex- 
act fixed point of the 2-loop FRG equation for A'(u), i.e. the 
derivative of (7) [95] using the value ( = e = 4 — d (here 
C — 4: which is the correct value for the present model). 
We can also check that (208) satisfies the general relation 



2A-(0+) _ g CT 



is) 



(209) 



In fact, we can also check some of the general relations dis- 
cussed in [62, 84] for the so-called higher Kolmogorov cumu- 
lants, defined there: 



G{X) := exp(A['u(w) —w — u{w') — w')] — 1 

= expl A] ) -1 



1 - yi - 4As^ 

2Sm 



-A) {w-w')+0{{w-w')'^) 

(210) 



for w > w' > wi, using formula (E6) for the Laplace trans- 
form of the first passage time probability with 6=1 and 
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D = .s,,,. This is related to the fact that the full distribution 
of avalanche sizes at the tree level (i.e. in mean field) in the 
field theory coincides with the distribution of return times of 
the 1-dimensional Brownian motion. 



B. Motion at finite velocity v > 

Let us now consider the case w{t) = vt, i.e. a particle 
puUed by a spring at constant velocity: 



dtu{t) = F{u{t)) + w?{vt - u{t)) . 



(211) 



For simplicity, we set 77 = l(it can be restored by setting t 
t/r] and v rjv). For v > Q and since it is an over-damped 
equation of motion (no overshoot), the instantaneous velocity 

vi := di ll is positive (possibly after a short transient), hence 
one can write Vt = y{u{t)), which satisfies: 



dtVt = F'{u{t))vt+m^{v-Vt) 



duy{u) 



m 



v(m) 



+ F'{u) 



(212) 
(213) 



Since for the Brownian-force landscape F'{u) is a white 
noise, the second equation is a standard Langevin equation, 
hence the probability to observe v at u, given that the velocity 
is Vq at uq satisfies the Fokker-Planck equation: 

duP = [ad,P + (9v^(v))P] (214) 
i?(v) = — m^f In V + m^v (215) 



with delta function initial condition. For u 
to the equilibrium measure: 



Qo it converges 



which yields the average velocity v*^ = v, as expected, 
and the connected expectation of the square of the velocity 

' = V:^. Note that the average velocity using P is 
= V + hence it does not even vanish as u 0+: 
this is because most (in fact, as m 0, all) of the u seg- 
ments belong to avalanches, yielding a finite average velocity 
if weighted by du, but that the fraction of time spent on them 
goes to zero, consistent with ^ 0+. 

ABBM also noted that the correlation function 

C{t-to) := j dvtdvo(vt-v){vo-v)Q{vt,t\vo,to)Qeq{vo) 

(221) 

satisfies the very simple equation (see also [85]): 

dtC = -m^C (222) 

obtained by multiplying (218) by {yt — v)(vq — v)Qeq{vo), 
integrating over vt = m and ?;o and using the fact that the 
current J = ad^vQ + m^(v — v)Q vanishes at v = and 
V ^ 00. After integration by parts one obtains (222). 

From this we can now obtain the renormalized disorder A 
at f > and discuss the crossover One first notes that from 
the definition (5), inserting w = vt and w' = vt' and taking 
two derivatives, one has 



C{t - t') = {vt - v)(vt' -v) = -v^m-^A"{v{t - t')) . 
Since from (222) and above, 

C{t) = wm-^e""'* , 



(223) 
(224) 



we obtain 



^-lg-E(v)/a 

vr{m'^v/a) 



v/a^—m v/cr 



One can also directly work with (212), rewriting it as a 
stochastic equation [77]: 



dvt = dF{t) +m^{v- vt)dt . 



(217) 



dF(ty = 2avtdt is a Brownian motion up to a time 
reparametrization, and the factor vt can be seen by writing 
Jdui^^f = Jdtvi\^)\ In Ito prescription this 
yields the Fokker-Planck equation for the probability Q{vt = 
V, t\vo, to) of velocity v: 

dtQ = d,[ad,{vQ) + m^{y-v)Q] (218) 
= [ad,{vQ) + \/{d,E{y))Q] (219) 



Hence the steady-state solution for f — > 00 is Qeq ~ v ^P^ 



the velocity factor originating from the change of variable 
from u to t. One has [77]: 



Qeq(v) = 



T{m'^v/a) 



-1+m vfa ^—m, y/a 



(220) 



A"(u;) = m^e" 

V 



-rn^w/v 



(216) Integrating twice we finally get: 



A(0) - A(w) =(jw- 



av 



(225) 



(226) 



a formula valid for any v > 0. The integration constant has 
been fixed by either of the two equivalent conditions: (i) no 
cusp at u! = 0; (ii) the large-w behavior is the same as in 
the bare model, and as in the statics, i.e. in the Umit v 0, 
given by (208). A non-zero velocity v > thus smoothens 
the cusp in a boundary layer of size w ^ v/m^, but the func- 
tion remains non-analytic: there is a subcusp, i.e. a non-zero 
A'"(0+); indeed one has at small w. 



A(0) - A(«;) = 



~2v' 



4 

ma 
'6^ 



+0{w'^) (227) 



This indicates continuity of u{w{t)) but jumps in its deriva- 
tive, the velocity. It remains to be understood whether this 
feature is more general or if it is tied to the long-range nature 
of the random force landscape. 

The distribution of avalanches times (in t) and sizes (in u) 
at 11 > can be extracted by studying the returns "near" the 
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origin of the process , i.e. the return to the origin inside the 
potential well -E(v). This is given by (212) or equivalently by 



(228) 



where dB{u) is a standard unit Brownian motion. Near 
V = we can first ignore the drift term — m^du. Define the 
change of variables y = 2au, then v(u) = r{y = 2au) = 



^i=i^i (y) is the norm of a d-dimensional BM in the vari- 
able y [92], which satisfies: 



dr ■ 



d- 1 



2r 
d=l + 



dy + dB{y) 



rn?v 



(229) 
(230) 



This yields the avalanche- size exponent r = 2 — rf/2, for 
d < 2, from the power law decay of first return probabiUtes of 
a Brownian near the origin, see appendix G 



_ 3 _ 
^ ~ 2~~2^ 



(231) 



below a critical velocity v < a/mP. This result was antici- 
pated in Ref. [90, 91]. Note that the definition of avalanches 
at > is not clear cut and requires a small velocity cutoff 
noted vq. In Appendix E 3 it is shown that: 
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(232) 
(233) 



for Vq ^ V. Since the drift is neglected (233) holds only for 
s ^ Sm = (r/m^ the large scale cutoff. The drift provides a 
large velocity cutoff v> = u/rn^ in (216), (220) and a large 



relaxation time cutoff 



n ^, with 



v>i>- 



For 



larger velocity v > a/m? the behavior changes qualitatively. 
It corresponds to d > 2, see appendix G, and the most proba- 
ble velocity in Qeq(v) is no longer near v = 0. 

More details and the solution including the drift term, are 
given in Appendix E 3 for the various regimes. In particular 
it is shown that one recovers the quasi-static size-distribution 
obtained above in the limit ?; = 0+ [96]. For v < a/m? 
the random walk v(m), in the continuum limit, comes back in- 
finitely often near the origin (i.e. near vq), hence the role of the 
drift term is mainly to cut off the rare large avalanches, very 
much like in the statics (see e.g. the discussion in [62, 84]). 
For V > G jrr? the random walk in velocity space is still 
certain to come back near the origin but only because of the 
drift. There are then two types of avalanches. In a fraction 
of them (computed in Appendix E 3) the instantaneous veloc- 
ity V does not reach v. these avalanches are still described by 
the model without the drift (first-return "time", conditioned to 
return) and lead to power law distributions. In the rest, the ve- 
locity reaches v and equilibrates in the well -E(v); the "time" 
between two returns at small vq <C can then be estimated as 
~ (vo)""* ""^"^ proportional to the inverse equilibrium prob- 
ability (either Peq or Qoq depending on whether one is inter- 
ested in avalanche size or duration). Typically there will be a 



bunch of small avalanches of the first kind separated by one of 
the second kind. Eventually at larger velocities returns to the 
origin vg <C w become very rare events and there is no real 
sense in which one can talk about avalanches. 



VII. DEPINNING AND EXTREME STATISTICS OF 
RECORDS 

A. Model without a mass: records without drifts 

For a particle pulled through a random-force landscape it is 
also possible to consider the problem without a parabola. The 
problem is easier to solve, but the correspondence with the 
FRG calculations is much less clear. Let us give here some 
elementary results. 

We now have to solve for the smallest root of 



(234) 



where / is the applied force. We study the case where the 
force is continually increased. The process then has 
jumps from Mj to Uj+i as the force crosses the values /j, 
which form an increasing sequence. These values are called 
the record values for the process F(u), and the Ui the record 
times. Statistics of records thus naturally occurs in the physics 
of depinning. The problem is to find the running maximum 
(i.e. the record) of an unbiased process, while in the case of a 
mass it had a drift. In the absence of a drift the only scale in 
the problem is the system size M. 



1. Uncorrelated forces 

Let us start with the discrete model of uncorrelated forces 
studied in Section III A, characterized by a force distribution 
Pf {F) for each site. A similar problem was studied in [72] 
(section IV-B). There the probability distribution of the full 
record value sequence (/i < /2, • • • < /„,...) was obtained 
for a semi-infinite line. It can be mapped onto a sum of in- 
dependent variables as follows: The sequence distribution can 
be obtained from 



*(/n)=E' 



(235) 



for any n > 1, where the are independent positive random 
variables, each with an exponential distribution P{a)da = 
e~°'da. The function: 



$(/) = - In 



dFPf{F) 



(236) 



describes the tail of the distribution, here the smallest /. For 

stretched exponential tails, as in class I, the growth is /„ ^ 
n^/'*^, while for power law tails, as in class II, the growth is 
exponential in n. 

Another set of results, remarkably universal, is known 
[73, 74] for the probability P{N\M) of the number of records 
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N, here equal to the number of jumps, for a system of size M 
(notations are inverted as compared to [73]). Then for an un- 
corrected sequence of Fi it was shown [86] that at large M 



TV = In M + ^VlnM , 



(237) 



where ^ is a univariate gaussian random variable [75, 76]. 
Hence the translationally-averaged avalanche size in absence 
of a mass should be 



1 ^ 
N ^ 



M 



M 
InM 



(238) 



Hence for large sizes M, the average number of jumps be- 
haves as [73]: 



N : 



M 



(240) 



There are also results for the jump sizes Si, named record 
ages li in Ref. [73]. The typical jump size is styp = 
M/N = V ttM /2 while the average maximal jump size is 
Sniax — 0.626508Af and the average minimal jump size is 
\/MA[73]. 



i.e. it is the typical avalanche size =. In the language of 
records the avalanche sizes Si = Ui+i — Ui are the time inter- 
vals between successive records, also called record ages. The 
translational average grows unboundedly with system size. 
Hence there are very few avalanches and they are almost as 
large as the system. Note that, at variance with the results on 
the sequence fi, this result is independent of the distribution 
Pf{F) for continuous distributions. 



2. Forces correlated as a random walk 



B. Model with a mass: records with a drift 

The usual problem of records with drifts [74] consists in 
studying the sequence: 



Yi — Xi 



(241) 



with c > where the Xi are symmetric random variables. 
One way to present the correspondence to the depinning 
model with a mass [97] is that 



In the case of a landscape obtained as a discrete-time ran- 
dom walk, Fi = Fi^i — iji where rji are uncorrected random 
variables drawn from the same symmetric continuous distri- 
bution P{—rj) — P{rj), it was recently obtained in Ref. [73] 
that 
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(239) 



Equivalently: 



M— ►oo " 



cj,{u) 




FIG. 18: 0(it) defined in (201) for /a random walk with D = 1.67x 
10^ (blue). The dashed (red) line denotes (f){u) = (walk with no 
drift). We show explicitly the graphical construction of u{w) for 
w = 3400, as well as a larger avalanche of size S — 4401. The 
system is discretized with stepsize 1 in ti direction. 



X, = ~F^, c^n? . (242) 

The set of (upper) records Yi^ , p — 1, . . . , N, i — 1, . . . , M, 
i.e. successive highest values, are the values rn^Wp at which 

ip+i occurs in the process 



a jump from Up — ip to Up+i 
u(w). 



1. Short-range correlations 

In the case of i.i.d. random variables with a drift it was 
shown that the total number of records N up to time M grows 
linearly as iV r{c)M with (normal) fluctuations which were 
characterized [74, 76]. However obtaining analytic results, 
even for r(c), for a general distribution was found difficult 
and some results were obtained only for special distributions 
P{F) [74, 76]. The function r(c) is related to the avalanche 
density 1/ (s) = r(c), which is finite in presence of a mass, 
and is computed here for small m? — c. We solved the prob- 
lem for arbitrary distributions P{F) and found universality in 
the small c-limit, with three classes. In addition we obtained 
the joint distribution P(w, s) of (i) the time s between one 
record and the next; (ii) the difference in value w with the 
previous record [98]. These results were given in in Section 
IV. 



2. Long-range correlations 

Let us now extend the discussion of Ref. [73] to records 
with drift, i.e. depinning with a mass. Again we consider the 
random walk Xi — + rji with i.i.d. random variables 

ri = —F of distribution Pf{ri). The alert reader will note that 
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Pf{ri), Pava("S) and P{s) below denote three different proba- 
bilities and functions. Pf {q) produces a correlated sequence 
Xn- Xo= For n > 1 we set: 

P(n) = Prob(X„ < Xq - cn) = Prob(y„ < Fq) 
Q{n) = Prob(Xi < Xq - ci,i = I, . . . ,n) (243) 
= Vroh{Yi<Yo,i = l,...,n) (244) 

The Sparre- Andersen theorem [87-89] states that 

Q{z) := f; Q{n)z^ = exp ( ^^") ' ^^^S) 



n=0 



n=l 



setting Q(0) = 1 by convention. We denote F{n) := 
(3(n — 1) — Q{n) the first passage probability that Yn crosses 
Yo between steps n — 1 and n. As in [73] the joint distribution 
of record ages (jump sizes) Si and number A'' of records is: 

P({s}, N\M) = F{si)F{s2) ■ ■ ■ F{sn-i)Q{sn)Sj^^ .^^^ 

' (246) 

While for c = 0, P{n) = 1/2 independent of n, leading 
to (239) and the very universal results of [73] quoted above, 
for c > the sequence P{n), hence Q{n), usually depends 
on the details of the distribution P{ri). Hence apart from the 
asymptotic behaviour at large n (hence M), one expects less 
universality. 

The following formula are still valid: The generating 
function for the probabiUty to have N records given M, 
P{N\M) = J2{s} ^({«}' ^1-^) can be written as 

OO 

P{N\M) = F{z)^-^Q{z) , (247) 

M=N-1 

where F{z) := En=i ^(")^" = 1 - (1 - z)Q{z). For 
instance the generating function for the average number of 
jumps is obtained by multiplying (247) by N, and summing 
over A'^: 
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M>0 



z^N{M) 



1 



{l-zYQ{z) 



(248) 



Similar results hold for higher moments. 
If one considers Pf {rj) with a finite second moment, the Xi 

are in the universality class of the Brownian motion and one 
should recover the results of Section VI, using that [89] 



P(n) 



~ n 



-l/2„-nS(c) 



Q{n) 



-3/2 -nS(c) 



(249) 



with a common function S{c) = 0{c?) at small c. 
For instance, if Pf{ri) is a univariate gaussian, P{n) = 

erf(cvW2)] ~ {2T;(?n)-^'^e-'''''/^ for large n. We 
will not study the Brownian case in detail, since it was al- 
ready discussed in Section VI, and we refer to Ref. [89] for a 
detailed asymptotic analysis (as well as a nice proof of (245)). 

Of course we expect that stable distributions play a special 
role. Here we detail one example of a Levy-type random-force 
landscape, for which the results for the records with drift are 



particularly simple, and present a nice generalization of Ref. 
[73], although they may not be as universal. Consider the 
Cauchy distribution, 



7r(?7^ + a?) 



(250) 



such that the distribution of X„ — is also Cauchy with 
parameter a na. Quite extraordinarily. 



P{n) 



! dx 



tt[x^ + {nay 
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7r[a;2 + (a)2 



--:p 



is independent of n, with < p = arctan(a/c)/7r < 1/2 for 
c> 0. Hence Q{z) = (1 - z)-p and F{z) = 1 - (1 - z^'P, 
and 



Q{n) = 



F{n) = {l-p) 



T{n + p) 
r(i + n)r(p) 

r{n + p-l) 



r(i + n)r(p) 



(251) 
(252) 



Using (248) one finds the average number of records (i.e. of 
jumps): 



— r{2 + M-p) M^-p 



r(i + M)r(2-p) 



r(2-p) 



(253) 



which grows as a power law of the size. Higher moments 
grow with the same scale: 



7V2 



-]V+ 2r(3 + M-2p) 



r(i + M)r(3 - 2p) 

Hence at large M the cormected fluctuations are 
2 1 



(254) 



r(3-2p) r(2-p)2 



]V^2(i-p) ^ (255) 



and in all cases the results of Ref. [73] are recovered for p = 
1/2, the case without drift. The full distribution takes a scaling 
form at large M: 



P{N\M) « MP-^gp{NMP-^) 



(256) 



Surmning (247) with this scaUng ansatz at large M, i.e. x :- 
— \nzKil — z small yields: 



Jo 



dM MP-^e-'^''gp{NMP-^) w x'Pe 



TT. 



with 31/2(1;) = e 

From (246) one sees that the distribution of avalanche sizes 

(i.e. record ages) is Pava(s) = F{s) for s = 1,2, For 

fixed p and large s it decays from (252) as a power law with 
r = 2 - p: 



Pava(s) = F(s) 



-(2-p) 



-r(p-i) • 



(257) 



This leads to a simple interpretation in terms of a directed 
random walk with traps of independent random release times 
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Si, distributed as P{s) - 8''^^+''^ and 1/2 < = 1 - p < 1. 
M is the total time t and TV the distance x traveled. As is 
well known, for fi < 1, x ^ and the distribution of z = 
t/x^/^ = Af/A^i/(i"P) is a Levy stable disti'ibution L^(z) 
with positive support, which is indeed the solution of (256), 

9piy) = fi-'y'-'/^'L^iy-^/^^). 

Although the strong universality of the symmetric case does 
not hold, we expect that all processes in the class of the 
Cauchy process remain critical even with a drift which has a 
power law distribution of avalanches given above, and a con- 
tinuously varying exponent. For stable processes intermediate 
between Cauchy and Gaussian, avalanches should be cut at a 
finite scale, which diverges with different exponents as c ^ 0. 
The situation of stable processes broader than the Cauchy dis- 
tribution remains open. 

VIII. GENERAL CONSIDERATIONS ABOUT AN 
Af-COMPONENT DISPLACEMENT FIELD 

Up to now, we have considered particles, and more gen- 
erally elastic objects and manifolds, whose position is a 1- 
component function. We now consider particles or elastic ob- 
jects embedded in higher dimensions iV > 1. For simplicity, 
we focus on a particle, but the considerations in this short sec- 
tion can be extended to d-dimensional elastic manifolds. In 
the next section we consider an application to a particle driven 
in a 2-dimensional random energy landscape. 

Consider two particles, which see the same random energy 
landscape, but which sit in different parabolas, labeled 1 and 
2. These parabola are chosen with the same curvature 
but their centers differ, and can have very different trajec- 
tories, which we call {■wi{t)} and {■W2{t)}. An interesting 
case is when the trajectories differ but the endpoints coincide 
""^1(^1,/) = W2{tij) as in the example of figure 19. In that 
example it is clear that 

\un,Aif)-UvjAts)\ ~ V2 A . (258) 

ky 




X 

FIG. 19: 2 particles dragged through a random energy landscape by 
parabolic potentials (springs) whose centers have identical starting 
and final positions but follow different paths {wi (t)} and {w2(t)}. 
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FIG. 20: 2 particles dragged through a random energy landscape by 
parabolic potentials (springs) whose centers follow parallel straight 
lines as described in the text 



This is because there should be a non-zero critical force, fc, 
and that each particle lags behind each parabola center Wi{t), 
up to fluctuations, roughly in minus the direction of drift. 
Thus the process it„, depends on the trajectory w{t), and to de- 
fine a single valued function u^, we have to restrict to a single 
well-defined trajectory w{t) in a quasi-static limit. Consider 
now figure 20. Both parabolas move with the same velocity 
V in x-direction. They are completely characterized by their 
position w — {x ~ xq + vt, y}, with w £ M^, y G R^^^. 
Especially note that without loss of generality, xq can be put 
to 0. Again we integrate the Langevin-equation (1), to define 
u{w). 

A more difficult question is whether u{w) depends on the 
initial condition u{to), and since we have set xq to 0, implic- 
itly on the starting time t^. One expects (see next Section how 
it occurs) that the dependence of ■u{w) on to disappears in the 
limit of to —00, and this is the limit we are interested in. 
This could be checked, similar to the exact sampling method, 
see e.g. [65], by starting at time to, and checking that at the 
time of interest t, all trajectories from all possible initial con- 
ditions have converged towards a single one. If not, one starts 
at an earlier time and checks again, repeating this proce- 
dure until all trajectories have converged. This defines a func- 
tion u{w), which is now independent of the initial time and 
conditions. In the next section it will be checked numerically 
that for a particle driven through a two dimensional bounded 
random energy landscape, all trajectories indeed converge, see 
figure 21. It is also found there that the so-called no crossing 
property (Middleton theorem [58]) does not hold for N = 2: 
although violations appear to be rare there are some instances 
of two trajectories splitting up. This results from a second par- 
ticle (more properly, a second trajectory of the same particle 
with a different initial condition) arriving at a later time on the 
same site: it then feels a stronger drive from the parabola and 
may jump forward and pass the first particle. An example is 
shown in figure 2 1 . 

Having given an unambiguous definition of u{w), we can 
calculate connected correlation functions of its moments. 
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FIG. 21: Trajectories of particles dragged from left to right by a parabolic well. Trajectories starting at different random initial positions all 
converge towards the same trajectory, with the same position at a given time (which is not visible on the plot). The minimum position of the 
parabola in y-direction is indicated by a straight (orange/grey) line. At the right, we show our coordinate system. The inset at the top is a 
blowup of part of the curve. It shows a deviation from the no crossing property which holds for A'^ = 1 (Middleton theorem [58]): Trajectories 
which are together can split up, even if later on they join again (see text). 



FIG. 22: Trajectories of 200 particles dragged from left to right. All particles start at (0,0). However, they sit in different potential wells (as 
described in the main text). As can be seen, they mostly move together on preferred trajectories, before separating again. The endpoints are 
joined by yellow lines, = 0.003. The particles in the two outer wells, as the center one, are marked in colors (red/blue, green). The well 
was moved for a total of 1000 steps. 



which again define A, (now a tensor), and higher cumulants. 
For this tensor the driving direction x will play a special role 
(for > 2 we expect isotropy in the other — 1 directions. 
The calculation is done for a particle in the next section. Other 
definitions of u{w) could of course be given. The simplest one 
is to pick a fixed but different driving direction. From statis- 
tical isotropy of the disorder the results should be the same 
up to the rotation. We defer the study of more complicated 
driving processes to future work. 



IX. A PARTICLE DRAGGED IN TWO DIMENSIONS: 

d = 0, TV = 2 

We now study particles dragged through a 2-dimensional 
random-energy landscape. 

The algorithm works as follows: We generate a random- 
energy landscape on a square lattice. A particle in addition 
sees a parabolic well. The total energy is the sum of both. 
We will mostly use a box distribution for the energy of a site, 
uniform in [0, 1]. Energies on different sites are uncorrected. 
We then update all particle positions: If a particle can move in 
a direction s.t. it will lower its potential energy, it will do so. 
If there are several such directions, it will choose the one with 



the lowest final energy. We allow moves to the eight nearest 
neighbors numbered from 1 to 8 (starting at the center 0); 
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If a move is possible, we perform it and then try other moves 
again, until the particle finally gets stuck. If several moves are 
possible, we take the one which results in the largest descent 
in energy, i.e. we go into the direction of the maximum force. 
Only then, we update the position of the parabola, by moving 
it from w to w + dw = w + vdt. We record the particle 
position as a function of time t. 

We first show numerically that there is a unique attractor 
trajectory (see figure 21). Start particles at random positions, 
but in the same parabolic well. Then move the parabola in 
a given direction (here always to the right, also denoted the 
x-direction). One sees that trajectories converge, and parti- 
cles will have the same position at a given time (not visible 
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on figure 21, which only shows the world-lines.) This conver- 
gence can be understood from the fact that if two "particles" 
(in fact these are the same particle but with different initial 
conditions) meet at a site at the same time, their future evolu- 
tion is identical. Hence the deep sites with low energies where 
the particle gets temporarily stuck act as sinks where the tra- 
jectories merge. Clearly, the particle needs to be trapped long 
enough for the process to be efficient. 

Our proper simulation is done with many particles (here 
200), each sitting in a parabola which are displaced by one 
unit (or in general by dy) to the top. This gives us data-points 
in the y-direction. In the x-direction (in which we move the 
parabolas), we use that after some time t the parabola has been 
displaced by a distance vt. We denote the minimum of the 
parabola wt^y and its zth component wl y as 

wt.,y ■■= {vt, y) , y := vt , wl y := y (259) 

The particle sitting in this parabola will have position u^_^ ^ , 
with components wj^^ . We then define A*^ {tv, y) as 



A^^- {{f - t)v, y) := «^ „ - wlo)«, , ~ ^l',y) 

(260) 

The connected symbol c indicates that we have subtracted the 
critical force. A*-' (a;, y) has the following symmetry proper- 
ties 

A'\x,y) = A"ix,-y) (261) 
A°\x,0)^ (262) 
A°i(0,y) = -Ai°(0,2/). (263) 

This is a consequence of the relabeling symmetry in (260) 
A^^{x,y) — A-''(— — y) and obvious covariance under 
the parity symmetry y —y. For smaller and smaller 
masses, there will be more and more data-points. Steps in 
the x-direction are necessarily discretized, of size dw = vdt. 
This poses an additional problem not present for = 1: 
there choosing a dw too large results in a loss in precision 
(since some smaller jumps may be overrun) but does not have 
dramatic consequences for large jumps, especially does not 
change the endpoint u{w), due to Middleton's theorem [58]. 
In contrast, for = 2, if the parabola is not moved adia- 
batically , the particles will see a strong force forward, and 
therefore be more likely to move forward, instead of sideward, 
thus embarking on a different trajectory. This may alter the 
whole trajectory over a much larger region. In practice, we 
decided to never move the parabola by more than one unit, 
before checking whether a move could be made. It may be a 
possible source for finite-size corrections. These will disap- 
pear if, and only if the critical force scales to zero for to ^ 0, 
since the energy gain for an elementary move is 




However, fc goes slowly towards 1, by which it is bounded. 
Unfortunately we find e.g. fdm^ = 0.01) = 0.464602, 
/^(to2 = 0.0001) = 0.785073. This might indicate that 




FIG. 23: Scaling collapse for Axx{x,0) :— 

m""''^'^''^ Axx{xm^^ ,0), with C,x ~ 1.595. The scaling col- 
lapse is perfect except for the two largest masses. 

the step-size we have used is still too large. We have not at- 
tempted to use a smaller step-size due to the enormous com- 
puting powers needed. We nevertheless believe that the results 
are valid for the following reason: fc measures the time aver- 
age of Uw — w, but we have to know the forward force exerted 
by the spring, when the particle arrives at the trap. Clearly, 
this must be much smaller, otherwise in a few steps the force 
would have increased by 1, which is sufficient to overcome 
any barrier for the box-distributed random energies, and the 
particle would not remain pinned for a long time. However 
we see diverging trapping times in the simulations, thus the 
argument using eq. (264) is not valid. 

We now present data for the force-force correlators in fig- 
ures 23 to 27, for masses ranging from = 0.1 to = 
10^^, descending in half-decades. A first and important qual- 
itative conclusion to be drawn is that all correlators not only 
depend on x, but also on y. This is in contradiction to the 
fixed-point structure used by Ertag and Kardar [63], whose 
Aij depends only on x but not on y. 

Our aim is to determine the scaling exponents C,x and C,y 
from the finite-mass scaling-ansatz, suggested by the FRG 
equations for this problem [66]: 

A]l{x,y) := 171-'^+'^'+'^' A'^ {xm-^- ,ym-^y) , (265) 

and supposing that Am A for to — > 0. We find that for 
A^^{x,0), Ayy{x,0), A^^(0,y), Ayy{0,y) and A„^(0,?;) 
separately such a scaling collapse is possible. There is no 
doubt that Cy = 1, with consistently rather small errors: the 
scatter from the different estimations is (y — 1.009 ± 0.015. 

However the results for are less consistent. We find 
different values, depending on which quantity we consider 
There is a clean data-collapse on figure 24 for Ayy with 
(^x — 2; however on figures 23 and 25 for A^x the best col- 
lapse is with Cxx = 1-6; finally on figure 27 for Axy, the best 
scaling collapse is for (x = 125, where however some of the 
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FIG. 24: Scaling collapse for Ayy{x,0) :— 

^-4+2Ca^^^(2,^C^ 0), with Co: = 2 and (y = 1. The de- 
scending lines indicate that no data has been collected for larger 
values, due to an insufficient choice of parameters. Good scaling 
collapse. 




FIG. 25: Scaling collapse for Axx{0,y) '■= 

m-'^+^'^^A^40,ym^y), with C:^ = 1.595 and Cy = 1- The 
scaling collapse is excellent for all but the three largest masses. 



data are noisy (but note that at least the data for the second- 
smallest mass, which are already very difficult to obtain, and 
for some of the larger masses show only little noise, s.t. noise 
does not seem to be an issue here.) Let us recall for com- 
parison that for iV = 1 a particle driven in a random energy 
landscape with a box distribution belongs to universality class 
III with a = 3, i.e. = 4/3, see sections IIIC 3 and HID 2. 
(For a box random /orce it is a = 1 and ( = 1). 

To conclude, we have shown that the unique attractor tra- 
jectory can be defined for the particle with N = 2, and that 
the disorder correlator Ay (x, y) can be measured. 




FIG. 26: Scaling collapse for Ayy{0,y) :— 

m"*^^'"" Ai:i(0, j/m^"), with C,y = 1. The collapse is good, 
except for the 3 largest masses. 




FIG. 27: Scaling collapse for AxyiOyU) := 

m-'^+'^^+'^yAxy{0,ym'^y), with = 1.25 and Cy = 1- 
The signal-to-noise ratio is rather big. In order to improve 
the statistics, we have used (after numerical verification) that 
A.y(0,y) = ~Ay40,y), to plot i[A,y(0,y) - Ay40,y)]. Fair 
scaling collapse, except for the 2 largest masses. 



X. CONCLUSION 

To conclude we have shown how the renormalized disor- 
der correlator A{w), central to the Functional RG theory of 
depinning, can be measured for a manifold of internal dimen- 
sion d driven by a spring in a iV = 1 random landscape. This 
correlator contains information about the stick slip motion of 
the interface. 

We have solved analytically the case d = of a particle in a 
short-range correlated pinning-force landscape, finding three 
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universality classes. In each case we have obtained the univer- 
sal fixed-point forms for A{w) for quasi-static driving, i.e. the 
depinning fixed point. In all cases it exhibits the famous cusp 
at w = 0+. We have also obtained the distribution of crit- 
ical forces, avalanche sizes and waiting times, and checked 
the general relations conjectured to hold between their mo- 
ments and the cusp A'(0+). While the exponent C can take 
various values depending on the class, the avalanche- size ex- 
ponent was found to be r = in all cases, which invalidates, 
at least at a naive level, the conjecture t = 2 — 2/{d + (). We 
also found that the distribution of avalanche sizes and waiting 
times are identical in the scaling limit. 

We have extended our results to a particle driven in a force 
landscapes with the correlations of a random walk. In the 
Brownian case, known as the ABBM model for interface mo- 
tion and Barkhausen noise, it is possible to solve for any driv- 
ing velocity and check the quasi-static limit. Remarkably this 
model has much in common with the mean field theory of 
avalanches and recent FRG results for avalanche distributions 
in d = 4. Since C = 4, the avalanche exponent for v = 
obeys, in that case, the conjecture t = 2 — 2/{d + (). 

In each case we have emphasized the cormections between 
the depinning problem for a particle and the extremal statistics 
of records, with and without drifts. 

These exact results in = provide interesting checks and 
interpretations of the Functional RG theory, and help us un- 
derstand what we should expect for manifolds. 

Finally, we started addressing the problem of depinning for 
iV > 1, mostly numerically; even for a particle it is quite non- 
trivial. We have checked numerically the consistency of the 
method, based on the ergodicity in presence of driving via a 
quadratic well. Numerous open problems remain. 
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APPENDIX A: DYNAMICAL ACTION AND OBSERVABLES 

In this Section we estabhsh the relation between the effec- 
tive action and observables; this is an extension to the dynam- 
ics of the proof given in [56, 57] for the statics. It. This rela- 
tion allows to measure the correlator of the dynamical FRG as 
explained in the text. Since the dynamical field theory is sig- 
nificantly more comphcated than the static replica field the- 
ory, the arguments presented here may be slightly less general 
and rely on further assumptions about the nature of the quasi- 
static limit and its ergodic properties, some remaining to be 



demonstrated, as e.g. extensions to > 1 components is less 
straightforward than in the statics [56, 57]. At a formal level 
however the arguments are rather similar. 

We use notations of the text and consider the equation of 
motion: 

T]dtu{x,t) = F^[u{t);w{t)]+^{x,t) (Al) 
Fx[u;w] = m'^{w{x) - u{x)) + Vlu{x) + F{x,u{x)) 

where w{t) is given, and ^{x,t) the thermal noise. We denote 
in general impUcitly = u{x, t;uo{x, to), F,^) the solution 
for given initial condition, disorder and thermal noise. 



1. Definition of functionals 

Let us first recall the definition of the useful functional of 
the dynamical field theory. We write the dynamical (MSR) 
action S in compact notation as: 

S[u, u\=u-g-'^ - u + u- A^°^ [u] - ■ B^^^ [u]-u. 

+0{u^) (A2) 

For any "vectors" u, v we denote u ■ v := j^^ UxtVxt (and 
additional index contraction for N > 1), A and B are respec- 
tively vector and matrix functionals. g^y can be an arbitrary 
(time independent) symmetric matrix but the usual choice is 
(in Fourier) g~^ = -{-m?. The functionals defining the bare 
action are 

A^^'^uU = Tldtu^t (A3) 
S^"' [u]xt,x't' = 2riT5.j;.j;'Stt' + Ao{uxt - Ux't')Sxx' (A4) 

which is the standard MSR action averaged over disorder. The 
statistical tilt symmetry (STS) of the bare action states that 

[u] and B^^^ [u] are invariant under the change u^t 
Uxt + (t^x- It implies the same symmetry for A[u\ and B[u]. 
The generating function of connected correlations is 

ex.\){W[w,'w]) := j 'D[u]'D[u\ exp — S[u,u\ 
+ u ■ g~^ -w + 11} ■ g~^ ■ u 



(A5) 



where the sources have been redefined, following [56, 57], in a 
convenient way for the following. It also admits an expansion 
similar to the action: 



W[w, w\= w ■ g ^ ■ w — w ■ A[w] + -w ■ B[w] ■ w 



(A6) 



where the STS implies that A[w] and B[w] are invariant under 
the change Wxt Wxt + 'f'x- Finally the effective action 
functional F is defined as usual as the Legendre transform of 
W: 



W[w, w] + T[u, u] = u- g 



w + w ■ 



u . 



(A7) 
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It admits the expansion: 

-1 



T[u, u] = u - g ^ ■ u + u- A[u\ — -u ■ B[u\ ■ u 



(A8) 



STS implies the same symmetry for A[u] and B[u] as for 
^(o)M andB(o)[u]. 



2. Relations to observables 

The functional W is directly related to observables, i.e. cor- 
relation functions, in the following way: Consider the average 
over solutions of the equation of motion, 



-S^[u,u\+w-g ^-u^t 



(eJ«t-9-^"')4,„„ = j V[u]V[u]c- 

= j D[M]D[ti]Me-'^["'"]+*-s~'-^+'^-s~'-" = e^[*"''^l (A9) 

Here S^, is the MSR action in presence of w, while S = Suj=q 
is the action defined above, in the absence of w. As usual, to 
take into account the initial conditions (if necessary) all time 
integrals start at to and an additional integral 2?[uo] -P[wo] fixes 
its probabihty at i = to. Expanding the above average, one 
finds: 



+^w-g ^ ■ {{uu)^- (u)^ {u)J-g ^-w 

(MO) 



Note that the (matrix) average (uu) — (u) (u) is not the 
connected thermal average but the connected double (disor- 
der+thermal) average. Note the index w which indicates that 
the above averages, e.g. {uxt)w, are averages w.r.t. Sw, i.e. in 
presence of the (given) driving Wxt- Comparing with (A6) we 
obtain: 



{Uxt)w = Wxt - gxyAyt [w] (All) 



{UxtUx't')w - {Uxt)w {Ux't')w = gxygx'y'Byt.y't'[w\ ,(A12) 

where summation (integration) over repeated indices is im- 
pUcit. Until now Wxt is arbitrary. For a uniform driving wt 
one has: 



{{wt - ut)) = m~ 



■'V 



(A13) 



{{Wt - Ut){wt> - Ut')) - {{Wt - Ut)) {{Wt> - Ut')) 

^ ' Byt^y't'[w] (A14) 



vv' 



These are the (spatially) local parts of the A and B function- 
als, and we expect: 

lim lim lim -ij / Ayt\w] = fc (A15) 



lim lim lim [ Byt,y't'\w] 

= L-''k{wt-Wt') 



More relations can be derived e.g. by considering the formal 
expansion (symbolically): 

+ 5w, w] 

= V -^{u..uu..u)^\g-Hw)^- {g-'w)^^ (A17) 
^-^ ni!n2! 

nin2 

Hence one has: 



{Ux't'Uxt)w = gxx'^tt' - {'^wA[w])x't',xt (A18) 

Let us also recall that the quadratic parts of W and F are (al- 
ways): 

Wqua.d[w, w]= w ■ g~'^ ■ R- g~'^ -w 

+lw-g-'-C-g-'-w (A19) 

rquad[M, U] = M • -R"^ • tt - • -D • U , (A20) 

where R and C are the exact response and correlation func- 
tions (in the absence of w) and C = R*DR. 

Finally it is useful to mention the terms without disorder: 

(A21) 
(A22) 
(A23) 
(A24) 
(A25) 
(A26) 

with: 

a-{l+g-a)=a (A27) 
a = a-{l+g- a)"^ (A28) 

3. Legendre transform 

Our aim is to relate the functionals A and B, which are 
observables as established above, to the functionals A and B 
associated to the effective action. 

Let us thus perform the Legendre transform. In this trans- 
formation one defines the functionals w[u,u] and w[u,u] 
which allow to compute T from W using (A7). One has 

w = g- — =u + g- A[u\ - g ■ B[u] ■ u + 0{u^) (A29) 



A[u] 


= a ■ u 


VuA[u] 


= a* 


g-^+a 


= R-' 


A[w] 


= d - w 


V^A[w] 


= a* 


g~^ - a 


= g-'-R-g 



Su 
ST 

w = g- — 

ou 



(A30) 



= u + g- WuA[u] -u - -g- V„u • B[u] ■ u + 0{u^) 
The inverse relations are useful as well: 

u = g- ^ =w~ g- A[w] + g ■ B\w\ ■ w + 0{w^) (A31) 

ow 



u = g- 



5W 
Sw 



(A16) 



(A32) 

3\ 



= w — g ■ VwA[w] ■ w + -g ■ ■ B[w] ■ w + 0{w ) 
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From (A29) and (A31) one deduces that 

w-u = g- A[u] + 0{u) = g ■ A[w] + 0{w) . (A33) 

This impUes that 

A[u] = iM| = A[u + g- A[u]] . (A34) 

which allows in principle to compute one functional from the 
other. One has the equivalent relation: 

A[w] =A[w-g- A[w]] . (A35) 

Moreover 

B[u] = -g-^ ■ dw/du , (A36) 

where u should be set to zero at the end. We have used a 
notation which makes the position of the indices clear. From 
(ABO): 

dw/dii = l + g- \/uA[u\ - g ■ V„u • B[u\ + 0{u^) (A37) 
Note the subtle difference with 

dw/du = l+g-A[u] Vm -g-B[u]-u^u +0{u'^) . (A38) 

This means that (at least at order 0{u) = 0{w), but eventu- 
ally even exact) 

g-'^ ■ dw/du = [g-'^ ■ dw/du]^ (A39) 
From (A31) follows 
B[u] = —g~^ ■ dw/du 

= -g~^d (u + g- A[w] - g ■ B[w] ■ w + 0{w'^j^ j du 

= d (-In + B\w\ ■ w + 0{iv^)j jdii 

= A\w\\J^ -g- B\u\^B\w\-dw/du , (A40) 
where it is implicit that u is set to zero at the end. It implies 

(l - [V^iM]* • 9)-B\u\ = B[w] {1+g- VuA[u])+0{u) . 

(A41) 

This can also be written as: 

(dw/dw)* • B[u] = B[w] ■ dw/du (A42) 

Equivalently: 

B[u] = (dw/dw,)* • B[w] ■ dw/du (A43) 

Finally, the relation between the B and B functional can be 
written as: 

B[u] = {l+g-VuA[u])'-B[w]-{l + g-VuA[u]) , (A44) 
where in this relation u and w are related via: 

w-u = g- A[u] = g ■ A[w] (A45) 
It can also be written equivalently as 

B[w] = (l-g- V^iH) ' ■ B[u] -{l-g- V„i[«;]) . 

(A46) 



4. Evaluation of the functionals 

Let us now evaluate the functionals Aa;t[{wj^f}] and 
Bxt,x't' [{"U^yt}] in various situations. 

Let us consider first T > equilibrium dynamics, i.e. a 
driving function Wyt which evolves infinitely slowly between 
Wyti = w-i{y) and Wyt^ = W2{y) such that the system always 
remains in equilibrium (i.e. we consider the limit f2 — i i oo 
at fixed wi — W2). From (All) and STS it is clear that 
A[w] = A[0] = in that limit. This impUes u = w in 
(A45) and also A[u\ = 0. It then implies that B[w] = B[w] 
and one recovers the results of Ref. [56, 57] for the statics 
using replicas. More precisely one expects in that hmit that 
Byt,y't'[w] = Byt.y't'[w] = d^^(^y)d^^(y,)R[wi,W2] where 

R = Ris the two-rephca functional of the statics. Hence it 
is a statement only about the infinitely separated time part of 
the B[w] functional and not about the smaller time separation 
part (which contains the renormalization of 77 and highly com- 
pUcated activated dynamics as described in [32]). 

Consider now T = and Wyt = w{t) = vt. From 
translational invariance Axt[w] (see e.g. (All)) can only be 
a time and space independent v dependent constant (assum- 
ing boundary conditions do not break translational invariance) 
which we choose to call f{m, v). Because of (A45) one must 
have the equaUty 

Axt[w] = Axt[u] = f{m, v) (A47) 

and w = w — m~^f{m, v) = vt — m~^f{m,v) in (A45). The 
difference with the equilibrium statics is that this constant is 
non-zero. This is allowed despite the STS symmetry because 
we are considering the T = limit first and the fact that w{t) 
depends on t cannot be ignored even for v = 0'^. In that limit 
one has fc{m, v — 0+) = fl = —fdm, v = 0^). Of course 
the fact that the constant depends uniquely on v assumes some 
ergodicity property, similarly if w{t) is a more comphcated 
adiabatic function there could be in general some history de- 
pendence. These issues have been discussed in Section VIII. 
For N = 1 we will rely on Middleton's theorem [58] which 
proves unicity of the solution. Note that at T > (A47) re- 
mains true with a /(m,f,T) such that f{m,,0,T) = in 
agreement with the discussion of the previous paragraph. 

Since the derivative of a constant is zero, using (A47) and 
(A44) we find that at T = with the choice w = vt one has: 

Byt,y't' [u = u{t) =vt- m~'^f{m, v)] 

= Byt,yt'[w]^A{v{t-t')). (A48) 

Denoting Byt,yt'[u] = A{u{t) - u{t')) gives 

A(u>) = A(w;) , (A49) 

a result on which is based our measurement of A(w) here and 
in [61]. 
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APPENDIX B: CALCULATION OF SOME INTEGRALS 

Here we compute the integrals in Eq. (84) of the main text. 
We need: 

/■OO 

/ dse-^lns = e-"lna-Ei(-a) (Bl) 

J a 

Ei(-a) = - / dte-yt (B2) 

J a 

This gives: 

A{W) 

pOO 

= (1 - e-^) / dae-«(i-<'"'") ln{a)e~^^'"' ln(ae-^) 

-(l-e-^) / dae-"(i-'^""')ln(a)Ei(-ae-^) 
Jo 

+6-^^(71 + y+^7£)-7l (B3) 
Consider the integral: 

- / dae-"(i-''""')ln(a)Ei(-ac-"') (B4) 
Jo 

POO POO 

= / da e-"^^-'^""') In a / die"*/* 

Jo Jae-W 



OO poo 



db / dae~''(i~'= ^ In 



1 ^0 

OO II POO 



a I die""' 



^ / dalnae-«(i--"^)— 
1 Jo 

/■°° d6 1 

A T (1 - e-^) + 6e-^ 

/■CSO 

X / dae-"(lna-ln((l-e-'^) + 6e-'^)) 



= -7b 



1 



X 



1 - e-^ 6(1 - e-w') 

2 ow2 , oi„2/ 



X 



TT- - 31^^ + 31n^(e^ - 1) + 6Li2 |^ — 
This yields 

1, 2 



1 



1 



A(a;) = y-^log^(e--l)-Li2(^^^) (B5) 
which can be rewritten as (85) in the main text. 



APPENDIX C: AVALANCHE PROCESS AND MARKOV 
CHAIN 

It is useful to recast the avalanche process for the discrete 
model of uncorrelated forces as a Markov chain, and define an 
algorithm for easy use in the numerics. 

Let us index jumps by n, they occur at positions u~ := w„ 
(integer), where the force is f„ (real). Note that the w- 
position of the jump us Wn = Un — m~^Fn. Given Fn) 



one finds the next jump (u„+i, F„+i) by the following algo- 
rithm: 

Un+l = Un + Sn , F^+l = /s„ (CI) 

Sn = min(p = 1,2,... such that fp<Fn + m^p) (C2) 

where /i,/2, •■ are a sequence of i.i.d. random variables of 
distribution Po(/)- The s„ (integers greater or equal to one) 
are the size of the avalanche and are determined at the same 
time. The variables (s„,i^„+i) form a Markov chain with 
conditional probabiUty P(s„, 



P{s, F'\F) = Po{F')9{F + m^s > F') JJ H{F + m^k) 

k=l 

(C3) 

s-l 

P{s\F) = {1- H{F + m^s)'[[H{F + m^k) (C4) 

pOC 

H{F) = dfPoif) 



k=l 



(C5) 



which is normalized X]is>i ^(^1-^) = 1 using that 1 — Hi + 
(1 - H2)Hi + (1 - H3)HiH2 + .. = 1 - H1H2H3... and the 
fact that the H). tend to zero as k increases. 

Starting from, for instance, Po(Fo), the distribution for P„ 
is given by: 

Pn{Fn) = P{Fn\Fn-l)P{Fn-l\Fn-2)..P{Fl\Fo)Po{Fo) 



P{F'\F)=Y,Pis,F'\F) 



(C6) 
(C7) 



S>1 



It converges to a stationary probability, noted P{F), which 
satisfies: 



^P(s,P'|P)P(P) = P(P') 



(C8) 



Once we find P(P) the joint distribution of (s„, P„_|_i, P„) is 
known: 

P(^^nj Pn+1 ; 

P„) = P(s„P„+i|P„)P(P„) (C9) 
Hence the avalanche size distribution is: 

P{s) = J dFP{s\F)P{F) (CIO) 
The sequence of waiting times is such that: 

W,, + l - Wn = Wn = Sn - m"^(P„+i - P„) (Cll) 

Hence the joint size and waiting time distribution is: 

P(W, s)=j dFdF'5{W -s + m-'^{F' - F)) 

X P{s,F'\F)P{F) (C12) 
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The problem is thus to determine the solution of (C8). One 
can formally write it as an infinite product: 



P(Fo) = f dFidF2.. f[ P{Fk-i\Fk) 
fc=i 



(C13) 



P{F'\F) = Po{F') J2 ^{F + m'^s > F') ]J H{F + m^k) 

(C14) 



S>1 



k=l 



which however also contains an infinite number of integra- 
tions. This method does not seem very practical (see however 
[76]) and in the text we obtain the result by another method. 



APPENDIX D: GENERAL RESULTS FOR 
UNCORRELATED DISORDER 



1. Proof that P{w) = P{s) 

Suppose, that the following variable transformation holds 
between j — w, and the corresponding 



f{j -w) = a^{j) 



(Dl) 



up to a constant shift and a rescaling. We recall that f(y) = 
exp(y) (Gumbel), f{y) = y", y > (Weibull), and f{y) = 
(-2/)"". 2/ < 0, (Frechet). 
Then 



P^iw',s>S) = j dyf"iy)e 



-fiy+W)-f{y+s)+f{y) 



(D2) 

where integration bounds depend on the class (real axis 
for Gumbel, positive axis for Weibul and negative axis for 
frechet). Using the relations in the text one finds that the joint 
waiting-time and avalanche size distribution is 



P{w >W,s>S) = 



pOO 

_ Jo 



/•OO 

/ <iyf"iy)e- 
Jo 



fiv) 



(D3) 

where the denominator is such that the distribution is prop- 
erly normalized. Also note that for some choice of variables, 
one has to be careful with the bounds of integration, see the 
Frechet classe, eq. (183). 

This formula shows that avalanche- size and waiting-time 
distribution are equal for all microscopic (uncorrelated) disor- 
der: P{S) = P{W). 



a. Distribution of local area a = ws 

The distribution of the local area a = ws defined in the text 
can be obtained as follows: 

OO OO 

P{a >A) = jdW jdS [d^d.Piw >W,s> S)] 0{WS - A) 



OO OO 

= jdW j dSP{w >W,s> S)dwdsO{WS - A) 







J(iW j dSP{w >W,s>S) 







X [SiA - WS) - WS5'{A - WS)] 



DC OC' 



-A— I dW I dSP{w >W,s> S)S{A - WS) 







(D4) 



We note a subtle point that when writing the last term as the 
derivative w.r.t. Aof-Aj P{w > W,s > S)6{A - WS), 
the explicit derivative of A cancels the first term. 

Inserting the integral representation (D3) for P{w > 
W,s> S) yields 

P{a >A) = 

OO OO 

Jdyf"{y) e-/(s'+^)-/(f+^/'^)+/(!') 







/ 

Jo 



dyf"{y)e 



-fiv) 



(D5) 



Note that boundary terms from the partial integration in (D4) 
vanish except for w = or s = 0, in which case A = 0. Thus 
we could possibly have a term ~ 6 {A). However we know 
that for w = we do not have a diverging probabiHty for an 
avalanche, and vice versa for w and s exchanged. The above 
result can also be written as 



P{a > A) 

_Ad 
dA 



OO OO 



dw dyf{y)e 



-fiv+e^)-fiv+Ae-n+f{v) 







/ 

Jo 



dyny)e 



(D6) 



For the Gumbel class, f{y) — e^, and this yields the formula 
given in the text. 



APPENDIX E: FIRST-PASSAGE TIMES AND AVALANCHE 
DISTRIBUTIONS 

In this Appendix we recall the basic method of the back- 
ward diffusion equation to compute distributions of first pas- 
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sage times on a half line and an interval, and then extend it 
to compute the first-passage-time distribution in the velocity- 
diffusion equation of the ABBM model. 



1. First-passage-time distribution on a Iialf line 

Let us call T{w'; w) the first-passage time at w' of a BM 
starting at position w < w' at time 0. Let us recall that the 
generating function 

G{w',p; w) = {e-Pn^''"^)) (El) 

satisfies the diffusion equation and boundary conditions 



DdlG + hd^G = pG 



(E2) 



G{w',p;w') = 1 , G{w', p = 0;w = -oo) = 

for 6 > 0. To see that, one introduces the diffusion kernel 
on X G [— oo,w'] in presence of an absorbing boundary at 
X = w', which satisfies: 

dtP = DdlP - bd^P = DdlP + bd^P (E3) 
P{x,t = Q\w,Q)=5{x-w) , P{w',t\w,0) = 

By definition of the exit time one has for p > 0: 

G{w',p;w) = - f dte-P'dt f P{x,t\w,0) 

Jo Jx<w' 

= l-p f Pix,p\w,0), (E4) 



with P the Laplace transform of P. The latter satisfies: 
pP{x,p\w,0)-S{x-w) = {Ddl+bd^)P{x,p\w,0) . (E5) 

Multiplying with —p and integrating over x from — oo to w' 
yields (E2). 

The solution of (E2) including the boundary conditions is 



G{w',p; w) = exp ( _ ^) ) . (E6) 



2D 



For 6 > it satisfies hnip^Q G{w,p;w') = 1. It is then 
inverted into the probability P{u; w' — w)du that T{w'] w) G 
[u, u + du[. 



P{u;W) 



W 



\/A-kD 



{bu - Wf 
4£)u 



e{u) 



(E7) 

where W = w' — w. Note that for negative drift 6 < one 
has: 



lim G{w,p; w') = 1 — q 
q := Prob(r(u''; w) = oo) = 1 — exp I — 



(E8) 



\b\{ w'-w) 



since in that case there is a finite probability q that the walk 
starting at w never hits x = w' . 



The Laplace transform of the probability to be at w at time 
t can be written as the probability to arrive there for the first 
time, and then repeatedly going with and against the drift, 
coming always back to w: 



LT,^^ ^'o__ ^-(t«-6f)V(4r>t) ^ 



Wo 



sJ^pD + 62 



b-\/4pD+b2 
-.e 2D ' 



2^^(«'--)i-nJ(.)n„-(.) 



(E9) 



where Hq (p) = Gb=|b|(M;o,_p; 0) is the return probability go- 
ing along the drift, and TV^^p) = Gb^^\b\(wQ,p;Q) going 
against the drift (with > a small cutoff which allows 
to cross in the microscopic model). Eq. (E9) expresses that 
the probability to be near a; = is a sum of n-th passage time 
events; the factor w\l2D = dt = d'W^/{2D) is the change of 
measure from time to space. 

Similarly consider the problem of the last passage time to 
of a Brownian at w = 0, with initial condition w = Oatt = 0. 
Its Laplace transform can again be expressed as a geometric 
series 

n,ast(p) = [1 - Ho (0)] , 

1 - nj(p)no {p) 

(ElO) 



Hence: 



niast(*o) 



{to) = 



2D 



{tl) 



12D^ 



(Ell) 
(E12) 



2. First-passage-time distribution on an interval 

Consider now a Brownian starting at w in an interval 
[wa,Wb]- Consider the functions Ga{w,p), Gb{w,p) and 
G{w,p) = Ga{w,p) + Gb{w,p) which satisfy the same dif- 
ferential equation (E2) but with boundary conditions: 

Ga{Wa,p) = 1 , Ga{Wb,p) = 

Gb{Wb,p) = l , Gb{Wa,p) = (E13) 

Then G{w,p) = (e"^^"''*^^^) is the generating function for 
the first exit time Tab{w) of the interval [wa, Wb] by a walker 
starting at w € [wa, Wb] at time zero. It satisfies (E2) with 
boundary conditions G(u;a,p) = 1 and G(u;f,,p) = 1. One 
finds: 



Gaiw,p) 



sinh(/i(t(;b - Wa)) 



Gb{w,p) = esc 



^ iwa-w) sinh(/x(w6 - w)) 
smh.{ij.{wb - Wa)) 



_ „^(u„-w) smh{fi{w - Wg)) 



(E14) 
(E15) 
(E16) 
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One checks that G{'w,p = 0) irrespective of the sign of b as 
the walk is certain to exit the interval. The function Ga{w,p) 
(resp. Gb{w,p)) is the same generating function restricted to 
walks exiting in Wa (resp. Wb), with normaUzations: 



Pa = Ga{w,p = 0) = 



exp(;^(wb -w))-l 
exp{-^{wb - Wa)) - 1 



Pb = Gb{w,p = 0) = l-pa 



(E17) 
(E18) 



for 6 > 0. 

One way to derive these results is to introduce the diffu- 
sion kernel P{x, t\w, 0) with absorbing boundary conditions 
at X — Wa and x = Wb- Its expression reads in Laplace: 

P{x,p\w,0) = 

^^(x-w) sinh{i^{x - Wa)) sinhji^jwh - wp)) 



= < 



g2o(x-TO) 



Z)/Usinh(;u(w6 - Wa)) 

ior Wa < X < w 
sinh(^('u;f, — a;)) sinh(/x(u' — Wa)) 

DiJ,smh{iJ.{wb - Wa)) 

forw < X < Wb 

(E19) 

It satisfies the diffusion equation with the proper boundary 
conditions. This is obtained as follows 



dte-P^dt / dxP{x,t\w,0) 



G{p,w) 



= D{dxP{x,p\w,0)\x=wa - dxP{x,p\w,0)\x=wJ 

(E20) 

where the first term is Ga{p,w) and the second Gb{p,w). 
They represent the flux from each boundary. 



3. Avalanche distribution for tlie ABBM model 

Let s = T(v', v) be the first passage "time" from v at "time" 
u = to v' < V at "time" s, the avalanche size (here de- 
fined with some velocity cutoff v'). The generating function 
G(v',p; v) = exp(— sp) satisfies the backward diffusion equa- 
tion: 



2 

ad^G + - m2)avG = pG 

G(v',p;v') = l G(v',p;+oo) =0 

a. Solution without the drift term 
From now on we denote: 



(E21) 
(E22) 



X = 



The solution, if one first drops the drift term is: 
G(v,p;v)=(-j 



(E23) 



(E24) 



One expects that it describes correctly small avalanches. 

There are two cases. For x < 1 one has G{p = 0) = 1, i.e. 
the velocity is certain to reach any fixed v' > 0. The leading 
behaviour at small v,v' is then for x < 1: 



G - 1 = -G(p/4cr)(i-^)/2 
C = 



(E25) 
(E26) 



Inverse Laplace transform from p to s yields a distribution 
l/s"^ with the value of r = (3 — .t)/2 given in the text. The 
small avalanche cutoff sq, necessary since is not nor- 
malizable at small s for r > 1, is provided by v,v' and its full 
form can in principle be obtained by inverse Laplace trans- 
form from p to s of (E24). Its order of magnitude is easy to 
read from (E24) as sq ~ v^/cr (see below a more precise esti- 
mate). 

One notes that the limit v' ^ can be taken for a; < 1 for 
any v. It is then easy to inverse Laplace transform (E24) in 
that Umit and to obtain the probabihty that if the velocity at 
u = is V, then the next stopping point u' is in the interval 
u' e [s, s -l- ds] as: 



G(0,5;v)rf5=-i^-(^re-W« 

r(^) s s 



so = vV(4a) 



M - (1 - •^)/2 



(E27) 

(E28) 



Next one can use the stationary distribution to find the proba- 
bility that choosing a u = the next stopping point is at s: 



(ivG(0,s;v)Peq(v) 



1 {iJ.)l [x)vs^^ a 



where to be consistent we have assumed s <^ Sm = cr/m'^. 

For X > 1, the probability of ever reaching v' > is G(p = 
0) = (^)^-i < 1, from (E24). At small v,v' (and of the same 
order) the following expansion holds for 1 < x < 3: 



G= (- 

V 



) V4cr 



(E29) 

This shows that, conditioned to returning near zero velocity 
the avalanches size distribution has again a tail fox 
So, with T = {1 + x)/2 > 1. 



b. Solution with the drift term 

Taking into account the drift term in (E21) reintroduces the 
large-size cutoff for avalanches at Sm ~ a/m^. 

It is also possible to solve the full problem, with the drift. 
One finds for the Laplace transform of the first-passage time: 



G(v',P;v) = (J)' 



(v-v') 



C/ l-vf-(l + ^^ 



),2-x,^^^^v 



),2-x, V^^v' 

(E30) 
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which reduces, for u = 0+ to the expression obtained for 
quasi-static avalanches: 



G(v',p;v) =e ^"'"^ 



(E31) 



from which, after inverse Laplace, (202) was obtained. 

One can now check that G{p = 0) = 1 for all x and v, v' > 
0; hence thanks to the drift the walk comes back infinitely 
often, as armounced in the text. 

Let us consider the leading behaviour at small v,v'. The 
expansion has the form: 



G = 



A{1 + 0(v^)) + Bv^-^ (1 + 0(v)) 
Ail + 0((v')2)) + B(v')i-'^(1 + 0(v')) 



where A and B are complicated functions of p. Here we focus 
on the case v < ajw?, then: 



G 



l = |(v-^_(v')-^) 



(E32) 



This yields 



G-\ 



r(l - + 

^ 2a ^„4+4p^i^ 



r(-i + ^) 
)) 



(E33) 



For p = m"^/ a one can check that this expression 

reproduces (E25) above hence in that case the distribution of 
small avalanches can be found neglecting the drift. From this 
expression Laplace inversion allows, in principle, to obtain the 
full avalanche distribution. We will not attempt to perform it 
but note that there is an additional pole structures for p = p„: 



Vn = 



n{n — x) 
which implies a decay 



n= 1,2,... , (E34) 



P{s) ~ e 



(E35) 



at large s:^ Sr, 



APPENDIX F: SOME 1-POINT OBSERVABLES FOR THE 
BROWNIAN FORCE LANDSCAPE 

We present here a few partial results for 1 -point observables 
for the Brownian force landscape, deferring a more complete 
study to the future. 

Note that since F{u) is an unbounded Brownian landscape, 
it has an infinite threshold force 2/(i = maxF(u)— minF(w). 
However the model studied here of a particle dragged by a 



parabolic well is well defined, and from it one defines an av- 
erage critical force: 



.fc{m) 



u{w)] 



m 



2-C 



(Fl) 



which diverges [99] as m — » for > 2. As discussed in the 
main text, for 7 = 1 one has C = 4. 

Let us now examine the 1 -point probability of the process 
u{w). Since it is a long-range correlated landscape, there is a 
subtlety linked to the choice of boundary conditions. 



1. Special boundary conditions 

If we first fix u{w = 0) = 0, the probability Pw{u) that 
u{iu) = M for M > is equal to the probability that the first 
passage time of (j){u') at w is u, starting at zero. Hence we 
find 



(F2) 



where P{u; w) is defined in (E7). In fact, it follows from the 
Markov property that, if we impose m(w = 0) = and leave 
the future unconstrained, the n-point probability is: 

Pwu-Wni'U'l, ■■,Un) = P{ui;Wi)P{U2 - Ui;W2 - Wi) 

X ... X P{Un - Un-i;Wn - Wn-l) 

(F3) 

for < wi < W2--- < Wn- Computing the moments of (F2) 
one finds: 



u{w) — w = 



2a 



{u{w) — wY = 2Dw = — jW 
mr 



(F4) 
(F5) 



hence one finds that the critical force fc{iTi), which is propor- 
tional to the average extension of the spring pulling the parti- 
cle is zerol On the other hand, the total area of the hysteresis 
loops per unit length cannot vanish - in fact from (F5) we can 
guess that it should grow as ^/w. Hence by contrast with the 
case of the uncorrelated force landscape, these two quantities 
cannot be equal. 

To understand this apparent paradox let us note that to in- 
sure u{w = 0) = one needs to impose rather strong condi- 
tions, e.g. F{u) = for all m < 0. Otherwise, there is a non- 
zero probabiUty that the BM has taken values ^(u) > w = 
in the past, i.e. for u < 0, which is in contradiction with 
u{w = 0) =0. If we now want to use the hysteresis loop ar- 
gument in a symmetric way, it would require a similar choice 
at some prescribed u{w = W) = 0. But then (F2) and (206) 
do not hold anymore (it holds for a Brownian unconstrained 
in the future). Since that procedure produces a non- stationary 
result, we do not discuss it further. 



2. A more generic situation 

First note that shifting F{u) by a constant leaves (200) un- 
changed, but also shifts w — u{w) while it does not change 
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the area of the hysteresis loop. Hence comparing fc{m) = 
m^lw — u{'w)] with the area of the hysteresis loop makes 
sense only for a statistically symmetric landscape. To get 
rid of this unimportant shift we can consider a distribution 
of forces symmetric around zero by setting F(u = 0) = 0. 
This is one way to fix the problematic zero mode of the Brow- 
nian landscape. This does not mean that Ui := u{w = 0) 
is necessarily at zero, ui is the position of the first crossing 
of 4){u) = w = by the BM. The probability distribution of 
Ml < is given by (see Appendix E): 



niast(wi) = , 

y 47rcr|ui 
This yields a critical force: 



^g-m^|«i|/(4a)^(_y^^ (F6) 



/c(m) = -rn^ui = 



2(7 



(F7) 



which makes more sense, i.e. it is positive and obeys the ex- 
pected scaling. 

However, to be a bit more general, what we have just com- 
puted is the conditional probability Pw^ {u\\F{ua) = fa) with 
the choice Ua = and fa = (which can be realized using 
shifts of the axis) andthe additional choice wi = Ua — fa/'m^ 
(hence wi =0). One wonders whether the critical force de- 
pends on that choice. 

To answer this question, we need to compute 
Pwi{ui\F{0) = 0) as a function of wi. There are two 
cases: If < then ui > and one has 

f«,i<o(wi|-F(0) = 0) 
= ^(-ui)LT;i,_„^Gi,=|i,|(|w;i|,p;0)niast(p) • (F8) 



This gives 
One finds 



/ 



d«iUiP^,<o(ui|F(0) = 0) = -^-^ (F9) 



If > there are two subcases, and the total reads: 

= 0{-Ui)LT~}^_^^Gb=-\b\{wi,p; 0)nia,st(p) 



x(l-Gb=_|i,|(w;i,p = 0;0)) 



(FIO) 



This gives: 



P^,>o{ui\F{0) ^ 0) 
\b\ 



0{-ui) 



V^nD\ui\ 



. ai \ "^1 -3/2 - 



(l-e-T^"'i) , (Fll) 



which can be checked to be correctly normalized. One gets: 

1° duiz.iP»,>o(wi|F(0) = 0) = - + ^) e-^-i 

duiuiP^,>oiui\F{0) = 0) = ^ (l - e-^'"^) 

Hence we find, setting |6| = 1 and D = a/m^, conditioned 
to F(0) = 0: 



'[wi — u{wi)] = 2a jm , wi < 



\w\ — u(w\)\ = 2{m wi + a/m )e 



2\ ~m wi/a 



(F12) 

, wi>0 



For the past {wi < 0) one recovers the previous result, while 
for the future it decreases to zero at large wi . 

We could try to generalize further by picking a Ua and av- 
eraging over fa with some distribution. Using the condition 
w = u— F{u) frr? one can check that: 

Pw^{u\\F{Ua) = fa) = P^i_„„+/„/m2(Mi - Ma|P(0) = 0) 

Hence we can use the previous calculation. For fixed w\, Ua 
and fa one has: 

j dUi{wi - Ui)Pyj^{Ui\F{Ua) = fa) 

= dui{wi - wi)P^j_„^+/^/„2(ui - Ua.\F{0) = 0) 

= y" du[iw[ - fain? - u[)P^,^{u[\F{0) = 0) 

{2cr/m2 - /„, for wi-Ua + fa/m? < 

2{w?w'^ + alm^)e-'^^<l'' - fa 
for w[ =wi-Ua + fa/rn^ > 

(F13) 

which can, in principle, be averaged on /„. It is not clear 
however at this stage which distribution to choose and how 
to relate these quantities to the area of the hysteresis loop. 
Further work is needed to clarify these issues. 



APPENDIX G: FIRST RETURN PROBABILITIES IN d 
DIMENSIONS 

In this appendix, we recall standard methods to derive the 
first return probability in d dimensions, using a regular hyper- 
cubic lattice. 

Consider a random walk on the integers Z. The probability 
to return to the origin after t steps (t = 2m even) is [93] 



Ut 



■ 

t/2 



(Gl) 



This is seen by noting that giving weight a for a step to the 
right, and h= 1 — a for a step to the left, the probabiUty to be 

at X after t steps is the same as taking i steps to the right and 
f — z to the left, with x = 2i — t and equal to px = a^b*~^ ( * ) , 
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from which the above result is obtained for a = 6 = 1/2, 
i = t/2. We note uq = 1. 

The probability to return to the origin for the first time after 
t steps is noted ft, and we note /o = 0. In order to be at the 
origin at time t, we must return for the first time no later than 
t, and can then make a new excursion. Therefore we have (for 
t > even) 



Ut = f2Ut-2 + fiUt-i + . . . + ftUQ 

Introducing the generating functions 

oo oo 
'^{x) := ^ U2mX^ , f{x) := ^ f2mX" 

m=0 m=0 

(G2) can be written as 

u{x) = 1 + f{x)u(x) . 
This gives for the probabilities (Gl) 
1 



f{x) = 1 - s/l^. 



thus the first return probabiUty in d = 1 at time t is 



ft 



At large times, 



U2n 



/-Km 



t- 1 



/2m 



(G2) 



(G3) 



(G4) 



(G5) 



(G6) 



(G7) 



We now want to calculate the same quantities in d dimensions, 
noted Ut and ff. uf is simply: 



= {utY 



As an example, in d = 2 



f=^x) = 1 



TT - ■ ■ 2K{x) 

with K the elliptic K function, and in d = 3: 



(G8) 



(G9) 



u{x) 



-\ 2 



(GIO) 



In general, an analytic solution for the coefficients is not pos- 
sible. We therefore give the asymptotic behavior for large t. 



i.e. small 1 — x, s.t. the series picks up contributions at large 
times. Noting s := — In x, we get for d < 2 



"2m''' 



m=0 



^"^7 — \d72 exp(-sm) 



^-rf/2gl-ir 1 



(Gil) 

Note that for d > 2, a UV cutoff is needed, which we discuss 
below. Using (G4) yields 



1 TT'^/^sl-^ 



'ix) r(i-l)- 

Transforming inverse Laplace gives 



(G12) 



J2i 



2-d d_ 



2m 



2„§-l . 



d-K 



2 -ms-'TTS"^ sm I — ) ford<2. (G13) 



Let us now consider d>2. Then 



u'^{x) ^ u'^{l) + T^-'^'^s^-^T 1 



(G14) 



u\l)^l + n-''/^(:(^ , d<2 (G15) 



«'^=3(1) 



= 1.3932 . . 



(G16) 



r(3/4)4 

where the approximation for u'^{l) valid for d near 2 was ob- 
tained summing using ~ 

gives 1.469 for 

d = 3). Then 



fix) = 1 



1 



u'^{x) 



1 - ^— + ^ TT-'^/^i-'T (l-- 

+0(s<^-2) . 



(G17) 



Note that 1 — is the probability that the walk never 

returns (which equals 0.2822 . . . in d = 3). Conditioned to 
returning, the probability that the first return occurs att = 2m 
steps decays at large-t as: 

u''{l)ftt2m « for d > 2 . (G18) 

We see that at d = 2 the exponent for the first-return probabil- 
ity changes direction, and that the amplitudes of both (G13) 
and (G18) go to 0. 
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